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ABSTRACT 

Existing  realization  procedures  for  a  fundamental  loop 
or  cut  set  matrix  are  reviewed,  compared,  and  classified 
broadly  on  the  basis  of  their  underlying  approach.   A  new 
combinatorial  synthesis  technique  is  presented  utilizing 
the  concepts  of  trunk  branches,  main  branches,  limbs,  and 
unique  connections  which  are  introduced.   This  procedure  is 
direct,  easy  to  apply  and  learn,  general,  and  yields  an  ex- 
pression for  the  number  of  physically  different  or  alternate 
realizations  which  are  possible.   A  general  computer  pro- 
gram for  realization  of  the  graphs  is  presented  and  illus- 
trated with  some  examples. 
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INTRODUCTION 

In  1847  Kirchhoff  presented  a  systematic  way  to  ana- 

0.9] 
lyze  complicated  electric  networks   .He  related  the  topo- 
logical structure,  circuit  elements,  voltages  and  currents 
in  the  most  general  form.   His  results,  except  for  the 
famous  "Kirchhoff s  Voltage  Law"  and  "Kirchhoff" s  Current 
Law",  were  almost  forgotten.   There  seemed  to  be  no  need 
for  further  investigations.   However,  during  the  last  ten 
years  engineers  started  to  investigate  properties  of  net- 
works which  previously  appeared  to  be  too  theoretical  and 
thus  of  no  practical  value.   The  reason  for  this  being 
that  many  modern  electric  network  problems  cannot  be  solved 
by  applying  conventional  analytical  methods  only.   Among 
these  problems  are:   The  generation  of  equivalent  networks 
with  a  minimum  number  of  circuit  components,  optimum  lay- 
out of  electric  networks  for  integrated  circuits,  and  the 
design  of  switching  networks  and  communication  nets. 

The  first  problem  that  was  attacked  by  several  inves- 
tigators was  the  realization  of  a  given  fundamental  loop  or 
cut  set  matrix.   These  matrices  play  an  important  role  in 
the  formulation  of  Kirchhoff 's  general  voltage  and  current 
laws  and  express  the  topological  properties  of  the  circuit. 
Since  this  theory  treats  the  electric  network  as  a  graph 
many  contributions  were  made  by  mathematicians.   Each  of 
the  existing  realization  procedures  in  the  literature 
suffers  from  some  of  the  following  limitations:   It  is 
either  too  abstract  and  cumbersome  to  be  of  practical  value, 
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it  involves  too  much  work  and  difficult  matrix  pattern  re- 
cognitions on  the  part  of  the  designer  to  be  of  value,  or 
it  fails  in  some  cases. 

In  this  thesis,  Chapter  2  presents  a  summary  of  the 
mathematical  representation  of  electric  networks  applying 
topological  formulations  as  currently  used  for  computer 
analysis.   In  Chapter  3  existing  procedures  for  the  realiza- 
tion of  a  fundamental  loop  or  cut  set  matrix  are  reviewed, 
compared,  and  classified  broadly  according  to  their  under- 
lying ideas  which  in  many  cases  were  determined  to  be  basi- 
cally similar,  a  fact  not  always  recognized  in  the  litera- 
ture.  In  Chapter  4  a  new  combinatorial  synthesis  is  pre- 
sented.  This  procedure  introduces  some  new  topological 
aspects  and  properties  of  a  graph.   Trunk  branches,  main 
branches,  limbs  and  uniquely  connected  branches  are  intro- 
duced.  The  realization  procedure  presented  is  not  only 
general  and  efficient,  but  also  provides  insight  into  the 
topological  properties  of  an  electric  network.   As  a  re- 
sult the  algorithms  developed  are  directly  applicable  to 
computer  programming.   A  detailed  general  realization  pro- 
gram (the  first  available  in  the  literature  to  the  author's 
knowledge)  is  presented  in  Chapter  5.   Moreover,  this  pro- 
cedure is  a  sufficient  condition  for  the  realizability  of  a 
loop  or  cut  set  matrix.   If  this  procedure  yields  a  graph, 
i.e.,  if  the  loop  or  cut  set  matrix  is  realizable,  the  to- 
tal number  of  physically  different  networks  corresponding 
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to  a  given  loop  or  cut  set  matrix  can  be  computed  as  pre- 
sented in  Chapter  6.   This  result  has  not  been  available 
previously. 
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MATHEMATICAL  REPRESENTATION  OF  ELECTRIC  NETWORKS 

2. 1   Introduction 

The  use  of  digital  computers  in  the  analysis  and 
design  of  electric  networks  has  not  only  provided  a  "super- 
sliderule"  for  numerical  calculations,  but  has  also  direct- 
ed the  attention  of  engineers  to  such  mathematical  areas  as 
linear  algebra,,  topology,  and  state  variable  theory.  Today 
classical  methods  such  as  Ohm's  law,  Kirchhoff's  laws,  and 
transform  theory  stated  in  terms  of  matrix  theory  and  topo- 
logy are  the  essential  tools  of  applied  electrical  engineer- 
ing, [l] 

In  this  chapter  it  will  be  shown   how  to  obtain  a 
graph  corresponding  to  an  electric  network.   The  most  im- 
portant notions  of  linear  graph  theory  will  then  be  defined 
and  the  concepts  of  fundamental  loop  and  cut  sets  introduced. 
Next  an  algebraic  description  of  a  graph  will  be  given  in 
terms  of  the  fundamental  loop  set  matrix,  B,  or  fundamental 
cut  set  matrix,  Q. 

Kirchhoff  introduced  many  of  these  ideas  in  his  origin- 
al work  on  networks  in  1847,  at  a  time  when  little  was  known 

[19] 
about  topology  „  During  the  last  fifty  years  his  results 

have  been  reexamined  and  generalized  by  men  like  Cauer, 

Foster,  and  Guillemin. 
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2 . 2  The  Electric  Network  as  a  Graph 

Given  an  electric  network  (Fig.  la)  made  up  of  lumped 
circuit  elements,  one  can  find  the  associated  linear  graph* 
(Fig.  lb)  by  the  following  rules: 

1.  Open  circuit  each  ideal  current  source. 

2.  Short  circuit  each  ideal  voltage  source. 

3.  Replace  the  passive  elements  by  edges  whose  end- 
points  are  called  nodes. 

4.  For  an  oriented  graph  (Fig.  lc) ,  redraw  the  graph 
maintaining  the  same  configuration  and  include 
the  directions  of  the  currents. 

2. 3  Description  of  a  Graph  in  Matrix  Form 

2.3.1  The  Incident  Matrix  A 

One  way  of  describing  the  graph  of  Fig.  1  is  to  list 
all  edges  and  nodes  and  to  specify  the  two  nodes  between 
which  each  edge  is  connected.   The  element  a. .  =1   if 

^  i: 

edge  i  is  incident   to   node  j.  When   the  graph  is  direc- 
ted, a.  .  =  -  1  depending  on  the  direction  associated  with 
edge  i.   By  definition,  a. .  =  +  1  if  the  ith  element  is 
directed  away  from  the  j th  node.   The  undirected  graph  of 
Fig.  1  is  listed  in  the  matrix  below. 

*    Since  most  terms  are  intuitively  familiar,  not  all  de- 
finitions are  stated  in  the  text,  but  are  listed  in 
Appendix  AI . 

**   See  Appendix  All  for  a  detailed  discussion  of  the  top- 
ological representation  of  sources. 
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edges 
12345678 
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1 

1 

0 

0 

1 

0 

0 

0 

b 

0 

1 

1 

0 

0 

0 

0 

1 
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0 

0 

0 

1 

0 

0 

1 

1 

d 

1 

0 

0 

0 

0 

1 

1 

0 

e 

0 

0 

1 

1 

1 

0 

0 

0 

nodes 


2.3.2    Trees    and   Co-trees 



A  graph  can  also  be  specified  by  its  edges  alone. 
Before  doing  so,  some  concepts  of  graph  theory  have  to  be 
introduced. 

A  graph  is  connected  if  there  is  at  least  one  path 
along  the  edges  between  any  two  nodes.   The  graph  of  Fig. 2c 
is  connected,  the  graphs  of  Figs.  2a  and  2b  are  not. 
The  unconnected  subgraphs  are  called  separate  parts.   A 
single  node  also  represents  a  separate  part  (Fig.  2b) .   An 
edge  is  removed  if  the  line  segment  between  the  two  nodes 
is  deleted,  the  two  nodes,  however,  remain  part  of  the  new 
graph.   Removing  edges  1  and  8  from  Fig.  lb  yields  the 
graph  in  Fig.  2c.   A  subgraph  of  a  graph  is  called  a  loop 
if 

(a)  the  subgraph  is  connected, 

(b)  exactly  two  edges  of  the  subgraph  are  incident 
with  each  node, 

(c)  exactly  two  nodes  of  the  subgraph  are  connected 
to  each  edge. 

Edges  (1,5,6)  of  Fig.  2a  form  a  loop 

Edges  (2,3,6,7,4,5)  of  Fig.  2c  do  not  form  a  loop 
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since  there  are  4  edges  incident  at  node  e. 
A  set  of  edges  of  a  connected  graph  is  called  a  cut 
set  (Figs.  3b  and  3c)  if 

(a)  removal  of  all  the  edges  of  the  cut  set  genera- 
tes two  separate  parts. 

(b)  Removal  of  all  but  one  element  of  the  cut  set 
leaves  the  remaining  graph  connected. 

We  can  now  introduce  the  concept  of  a  tree.  A  tree 
of  a  connected  graph  is  a  connected  subgraph  which  con- 
tains all  the  nodes  of  the  graph,  but  does  not  have  any 
loops. 

The  edges  of  the  graph  which  are  contained  in  the  tree 
are  called  tree  branches  or  just  branches. 

The  edges  of  the  graph  which  are  not  contained  in  the  tree 
form  the  co-tree.   They  are  called  chords  or  links . 
The  same  graph  (Fig.  lb)  can  have  different  trees  (Figs. 
4a-c) .   The  graph  in  Fig.  4d  does  not  represent  a  tree. 
Once  we  have  specified  a  tree  (T)  for  a  connected  graph  of 
n  nodes  and  e  edges,  then  following  properties  hold: 

(a)  There  is  a  unique  path  along  the  tree  between 
any  pair  of  nodes. 

(b)  There  are  b  =  (n-1)  tree  branches  and  l  =  e-(n-l) 
links,  where  n  is  the  number  of  nodes  and  e  the 
number  of  edges. 

(c)  Every  link  of  T  and  the  tree  path  between  its 
nodes  form  a  unique  loop  (fundamental  loop) . 
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(d)   Every  cut  of  the  graph,  which  besides  links  cuts 
only  one  tree  branch,  forms  a  unique  cut  set  of 
the  graph  (fundamental  cut  set) , 
2,3.3  The  Fundamental  Loop  Set  Matrix  B  and  the  Fundamental 
Cut  Set  Matrix  Q 
To  find  B  and  Q  for  the  graph  in  Fig.  lb  we  first 
choose  a  tree,-  say  that  of  Fig.  4b.   The  columns  of  the 
fundamental  loop  set  matrix,  B,  correspond  to  the  edges  of 
the  graph,-  and  the  rows  to  the  fundamental  loops  (which  are 
determined  by  the  links) . 

The  reference  directions  of  each  fundamental  loop 
shall  be  that  of  the  defining  link  as  indicated  in  Fig. 5  . 
If  edge  k  is  not  in  loop  i,  b.,  =  0   If  the  direction  of 
edge  k  is  opposite  to  the  reference  direction  of  loop  i, 


then  b.n  =  -1 
ik 


B   = 


For  Fig.  5a 
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t- 
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0 
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0 

0 

-1 

1 

1 

3 

0 

0 

0 

1 

0 

-1 

1 

0 

4 

I    loops   { 1.1) 


links        branches 
B  can  always  be  partitioned  as 


B  =  L-i 


U„   F 


(1  2) 


Where  U„  is  the  I   x  i  identity  matrix  and  F  is  a  b  x  i 
matrix.   The  columns  of  the  fundamental  cut  set  matrix,  Q 
correspond  to  the  edges  of  the  graph  (as  was  the  case  in 
the  B  matrix) ;   the  rows,  however,  correspond  to  the  fun- 
damental cut  sets  (determined  by  the  tree  branches) .   The 
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(a) 


(b) 


(c) 


Figure  1 

An  electric  network  (a)  and  the  corresponding  undirected  (b) 
and  directed  graphs  (c) . 
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(a) 


(b) 

Figure  2 


(c) 


Separate  parts  (a,b)  of  a  graph,  subgraph  (c)  which  is  not 
a  loop. 


(a) 


(c) 


(b) 

Figure  3 
Cut  sets,  #1  in  (b)  and  #2  in  (c) ,  of  a  graph  (a). 
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(c) 


(d) 


Figure  4 


Trees  (a,b,c)  of  graph  in  Fig.l,  subgraph  (d)  which  is 
not  a  tree. 


(a)  (b) 

Figure  5 

Fundamental  loops  (a)  and  cut  sets  (b)  of  a  graph. 
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reference  direction  of  each  fundamental  cut  set  shall  be 
that  of  the  defining  tree  branch  as  indicated  by  the  short 
arrows  in  Fig. 5b.  If  edge  k  is  not  in  the  cut  set  i,q.,  =  0, 
if  the  direction  of  edge  k  is  opposite  to  the  reference 
direction  of  cut  set  i,  q.,  =  -1.   For  Fig.  5b 
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Q 


links         branches 
Q  can  always  be  partitioned  as 


b  cut  set   (1.3) 


Q  =  L 


=   E  ,  U. 


bJ 


(1.4) 


Following  relation  holds  (see  Section  2.4.2.  for  the  proof). 

T 
E  =  -  F 

2 .4   General  Network  Analysis 

In  classical  network  analysis  two  methods  serve  to 
analyze  an  electric  network.   For  loop  analysis  the  loop 
currents  are  chosen  as  variables  and  Kirchhoff 's  Voltage 
Law  (KVL)  is  applied  to  each  loop.   For  nodal  analysis  the 
node-to-datum  voltages  are  chosen  as  variables  and  Kirch- 
hoff s  Current  Law  (KCL)  is  applied  to  the  n  nodes. 

These  classical  methods  shall  be  generalized  and  a 
third  method,  the  state  variable  method,  be  introduced. 
To  simplify  the  presentation  we  will  assume  that  the  imped- 
ance of  each  edge  is  resistive. 
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2.4.1   Loop  Analysis 

The  response  of  a  network  is  completely  known  if  the 
voltages  across  and  the  currents  in  all  its  e  edges  are 
known.   If,  however,  the  impedance  of  the  edges  are  known, 
then  either  the  set  of  e  edge  voltages  or  e  edge  currents 
suffices,  since  these  two  sets  are  related  by  Ohm's  Law. 

Once  a  tree  is  chosen,  fewer  than  e  variables  suf- 
fice to  characterize  the  network.   It  has  been  shown  that 
the  removal  or  opening  of  all  links  does  not  leave  any 
loops.   Expressed  in  circuit  terminology,  setting  the 
link  currents  equal  to  zero  forces  all  the  tree  branch 
currents  to  be  zero  (KCL) .   Hence  the  set  of  I    link  cur- 
rent (j)  can  be  represented  as  a  linear  combination  of  the 
link  currents (i). 

It  can  be  seen  by  inspection  that 

2  =  BTi  (1.5) 

where  the  components  of  the  vectors  j  and  i_  correspond  to 
edge  and  link  currents  respectively. 

Applying  KVL  to  each  loop  the  following  equation  is  obtain- 
ed 

B v  =  0  (1.6) 

where  the  components  of  v  correspond  to  edge  voltages. 
Representing  a  typical  edge  as  in  Fig.  6  and  applying  Ohm's 
Law  equation  (1.7)  is  obtained. 

v  =  Rj_  +  vg  -  Rj_s  (1.7) 
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where  R  is  an  (e  x  e)  diagonal  matrix,  r..  representing 
the  resistance  of  the  ith  edge,  and  r.  .  =  0  for  i  ^  i . 
j,    j  /  v,  and  v   are  (e  x  1)  column  vectors,  whose  ith 
component  represents  current,  current  source,  voltage, 
and  voltage  source,  respectively,  of  the  ith  edge. 

Combining  equations  (1.5),  (1.6),  and  (1.7)  yields 

BRBfc  i  =  -  Bv   +  BRj 

or 

Z^i  =  es  (1.8) 


where 


Z£      =    BRBfc  (1.9) 


is  called  the  loop  impedance  matrix,  and  is  of  order  i 
and 

es  =  -  Bvs  +  BRj_s  (1.10) 

is  called  the  loop  voltage  vector  {I   x  1) . 

2.4.2   Cut  Set  Analysis 

It  has  been  shown  that  a  tree  connects  all  the  nodes. 
If  the  tree  branch  voltages  are  forced  to  zero  (short  cir- 
cuiting the  tree  branches) ,  then  all  nodes  coalesce,  forc- 
ing all  edge  voltages  to  zero.  Thus  the  set  of  b  =  (n-1) 
tree  branch  voltages  is  the  only  independent  set  of  edge 
voltages,  i.e.  each  edge  voltage  (v)  can  be  represented  as 
a  linear  combination  of  the  tree  branch  voltages (e). 

It  can  be  seen  by  inspection  that 

v  =  Qte  (1.11) 
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where  the  components  of  the  vectors  v  and  e_  correspond  to 
edge  and  tree  branch  voltages  respectively. 
Applying  KCL  to  each  cut  set, 

Qj_  =  0  .  (1-12) 

Using  Ohm's  Law  for  each  edge, 

j  =  Gv  +  j   -  Gv  (1.13) 

—    -"  s     — s 

where  G  is  an  (e  x  e)  diagonal  matrix,  g. .  representing  the 

conductance  of  the  ith  edge,  g. .  =  0  for  i  ^  j.   j,j   and 

1D  —  — s 

v,v   represent  the  currents,  voltages  and  sources  of  each 

edge  as  in  equation  (1.6). 

Combining  equations  (1.11),  (1.12),  (1.13)  yields 

QGQfce  =  QGvs  -  Qj_g 


or 


where 


Yqe=is  (1.15) 


Y   =  QGQfc  (1.16) 

q 

is  called  the  cut  set  admittance  matrix  and 

is  =  QGvs  -  Qjg  (1.17) 

is  called  the  cut  set  current  source  vector. 

Properties  of  the  B  and  Q  matrices  which  are  important  for 

the  realization  procedures  described  later  in  this  thesis 

follow. 


Theorem  1 : 


QBfc  =  0  (1.18) 
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Proof:   From  equations  (1.5)  and  (1.10)  we  have 


Thus 


t 
Bv  =  0   and   v  =  Q  e 


BQfce  =  0 


Setting  any  e.  =  1  and  all  the  other  components  of  e_ 
equal  to  zero  requires  that  each  element  of  BQ   equals 
zero. 

Using  equations  (1.2)  and  (1.4)  the  following  rela- 
tion holds : 


Theorem 

2 

: 

E  =  -  FL 

Proof: 

QBfc  = 

,.E  > 

1 

E 

+  Ffc  = 

3 

E  !ubJ  &-J  aa  <ia9> 

F 


Thus  the  fundamental  loop  and  cut  set  matrices  can  be  ob- 
tained from  one  another  by  the  relation 

_U  1  Fj   and   Q  =  [-F*1  |  UbJ  (1.20) 


B  = 


2.4.3   State  Variable  Analysis 

The  description  of  electric  networks  by  state  variables 
is  included  here  because  its  formulation  rests  heavily  on 
the  concept  of  trees,  co-trees,  and  of  network  graphs. 
The  advantages  of  state-variable  analysis  are  [4]: 

1.  It  represents  a  unified  approach  to  networks 
having  non-linear  and  time-varying  elements. 

2.  It  provides  insight  in  the  behavior  of  networks 
in  such  areas  as  sensitivity  and  stability. 
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3.  It  is  better  suited  for  digital  computation  since 
an  nth  order  differential  equation  is  represented 
by  a  set  of  n  first  order  differential  equations. 

Definition:  [5]   The  state  of  a  system  is  the  minimum  set 

of  variables,  the  state  variables,  which 
contain  sufficient  information  about  the 
past  history  of  the  system  to  permit  the 
computation  of  all  future  states,  given 
that  the  future  inputs  and  the  state 
equations  of  the  system  are  known. 

This  definition  does  not  attach  any  physical  significance 
to  a  state  variable,  and  it  is  up  to  the  engineer  to  make 
a  proper  choice.    One  choice  of  state  variables  for  an 
RLC  network  consists  of  capacitance  voltages  and  inductance 
currents  (assuming  there  are  no  loops  consisting  of  capaci- 
tances only  and  no  cut  sets,  which  are  made  up  completely 
of  inductances)  [6]: 

This  choice  seems  quite  legitimate,  since  independent 
inductors  and  capacitors  are  the  only  energy  storing  ele- 
ments.  An  example  shall  illustrate  how  the  state  equations 
are  found.   A  linear,  time-invariant  network  (Fig.  7a)  is 
chosen,  which  does  not  contain  C-loops,  nor  L-cut  sets. 
The  method  can  be  modified  to  be  applicable  to  these  cases 
[7]. 

Choose  a  tree  (Fig.  7b)  that  includes  all  capacitive  edges 
and  as  many  resistive  edges  as  necessary.   This  tree  is 
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called  a  proper  tree  [6].   The  links  then  will  be  made 
of  inductive  edges  and  the  remaining  resistive  edges. 
Using  KVL,  Bv  =  0_,  it  is  convenient  to  partition 

v 


up 


-  Lvni  vx i v_ i v_J  ,  where  the  subvector 
— R  i  — L  i  — C  i  G 


s  vR,  vL#  v^,  and 


v_  correspond  to  the  voltages  of  resistive  and  inductive 
links  and  capacitive  and  resistive  tree  branches,  respect- 
ively.  Thus 

'"v. 


Bv  =  0 


or 


u*:fj 


^G 


=  0 


For  the  example  given: 


^R 


B   = 


v 


Rl 


v 


V 


SI 


R2 


1 

2 

3 

4 

5 

6 

7 

8 

1 

0 

0   0  , 

1 

1 

0    0 

1 

0 

1 

0   0  ' 

1 

1 

1    1 

2 

0 

0 

1   0 

0 

-1 

-1   -1 

3 

0 

0 

0    1   ! 

i 

0 

-1 

-1    0 

4 

VL3 

VC5 

VR7 

^L 

VL4 

^C  = 

VC6 

^G  = 

VR8 

_      

Using  the  same  partitioning  we  write  KCL  as 


-Ffc  1  U, 


2c 

2g 


=  0 
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Figure  6 


Typical  branch  of  an  electric  network 
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Figure    7 


Proper  tree  for  state  variable  representation 
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For  the  example  given 


2r  = 


Rl 


R2 


and  j  ,  j  ,  j   having  the  same  form  as 
— J_i  — (_   — (j 

vt  >    Y_,  Y_  substituting  j  for  v. 


To  eliminate  the  unwanted  variables  (v  ,  v_  ,  v_,  j_,  j  , 

R    Li    G   — R   — C 


i  )  we  apply  Ohm's  Law: 


XR  =  RlR  +  % 


2g  =  G^G  +  ±g 


^l  =  L  dt  1L  +  * 


2c  =  c  dt  ^c  +  ^c 


Where  e_  and  _i  stand  for  independent  sources 


R  = 


Rl 

0 

L   = 

L3 

0 

G   = 

G7 

0 

c  = 

C5 

0 

0 

\ 

0 

Li 

0 

\ 

0 

°t 

Combining  KVL,  KCL  and  Ohm's  Law,  we  obtain  the  state  vari- 
able equations  in  normal  form 


C 

0 

^C 

d 

dt 

__0 

L 

A 

All   A12 


A      A 

21    22 


pc 

ic 

+    B 

% 

(1.21) 


where  A. .  and  B  are  expressible  in  closed  form  in  terms  of 
the  matrices  F,  R,  L,  G,  and  C  [4]. 

The  solution  of  the  normal  form  (1.21)  can  be  found  in 
a  number  of  ways  using  a  digital  computer  [8],   Any  set  (y) 
of  currents  or  voltages  of  the  network  can  then  be  expressed 
in  terms  of  the  state  variables  (x)  and  the  sources  (u)  as: 
Y_  =  Dx  +  Eu  (1.22) 
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REVIEW  OF  EXISTING  REALIZATION  PROCEDURES  FOR  THE 
FUNDAMENTAL  LOOP  AND  CUT  SET  MATRIX 

3. 1   The  Realization  Problem 

It  has  been  shown  that  the  fundamental  loop  set  matrix 
B  and  cut  set  matrix  Q  are  essential  in  the  general  analy- 
sis of  electric  networks.  Moreover,  switching  circuits  [15] 
and  communication  nets  [3]  can  also  be  described  convenient- 
ly by  these  matrices.   A  next  logical  step  is  to  use  the 
loop  or  cut  set  matrix  in  a  synthesis  procedure  d 5]  and  [3]. 

It  was  relatively  easy  to  find  the  B  and  Q  matrix  for 
a  given  network.   To  realize  a  graph  corresponding  to  B  or 
Qf however,  is  not  a  trivial  problem  because  there  may 
exist  no  graph  at  all,  there  may  exist  exactly  one  graph, 
or  there  may  exist  many  graphs  having  the  same  B  or  Q  ma- 
trix.  An  exhaustive  search  procedure  forbids  itself.  Since 
1958  [16]  several  systematic  synthesis  procedures  have  been 
developed,  which  can  be  classified  broadly  according  to 
their  underlying  ideas.   In  this  chapter  the  same  synthe- 
sis problem  is  solved  using  a  method  typical  of  each  clas- 
sification, so  as  to  compare  the  procedures  and  the  work 
required.   The  important  theorems  used  in  these  procedures 
are  illustrated.   Their  exact  proofs  can  be  found  in  the 
references  as  noted. 
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3 . 2   The  Linear  Tree  Port  Structure  Approach 

Using  the  properties  of  the  admittance  matrix  Y  of 
order  n  of  a  linear  tree  port  structure  ,  Biorci  [9]  and 
Kim  [lo]  find  the  corresponding  graph  having  (n  +  1)  nodes 
and  n  ports.   Since  Y  =  QGQ   (1.16),  this  procedure  can  be 
used  to  realize  a  fundamental  cut  set  matrix  by  setting  all 
edge  conductances  equal  to  one.   In  Fig.  8,  using  edges  1, 
2,  and  3  as  tree  branches  for  ports  1,  2,  and  3,  the  corre- 
sponding fundamental  cut  set  matrix  is 


Q  = 


11 

12 

13 

22 

23 

33 

1 

1 

1 

0 

0 

0 

0 

1 

1 

1 

1 

0 

0 

0 

1 

0 

1 

1 

Forming  Y  =  QGQ   the  elements  of  row  1  of  Y  are  y, ,  =  g.   + 

g12  +  gir  y12  =  g13  +  g12.  y13  =  g13. 

Assuming  positive  edge  conductance,  we  observe,  that 
the  value  of  the  elements  of  Y  taper  off  to  the  right. 
Thus 


Y 


11  S  y12 


y 


In 


y22  * 


y 


2n 


Similarly 


y 


nn 


y 


n-l,  n 

y  -.  a  y  ~> 

^n-l,n  Jn-2,n 


In 
n-l,n 


*   Tree  branches  from  a  linear  path  and  the  terminals  of 
each  port  are  identified  by  the  end  nodes  of  a  tree 
branch. 
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These  observations  can  be  generalized  [3]  into  the 
following : 

Theorem  3 :   The  admittance  matrix  Y  of  a  resistive  n-port 
(with  a  linear  tree  port  structure)  in  which 
the  polarity  of  the  port  voltage  and  ordering 
of  each  port  is  aligned  in  one  direction 
(Fig.  8) ,  is  a  uniformly  tapered  matrix. 

A  matrix  is  uniformly  tapered  if 

1.  Each  element  is  greater  than  or  equal  to  zero 

y\.  £  0,  for  all  i,  j, 

2.  The  magnitudes  of  the  elements  of  each  row  of  Y 
taper  off 

(a)  from  the  main  diagonal  to  the  right  hand  side 

y .  .  ^  y .       for  j  ^  i 
J ij        J 1, 3+1 

(b)  from  the  main  diagonal  to  the  top,  the 
diagonal  element  being  the  largest  one. 

y . .  ^  y .  ,  . /  for  1  <  i  s  j 
i]    ^1-1,3 

3.  y.  .  +  y.  -,     ■    -,    ^  y ■  n  .  +  y.  ■  ,  -,    for  all  i  £  j 

Kim's  method  works  best  for  a  linear  tree,  Biorci's  method 
works  more  efficient  if  the  tree  does  not  form  a  linear 
path.   Kim's  method  is  used  to  solve  an  example.   His  real- 
ization procedure  consists  of  four  steps: 

1.  Form  Y  =  QQt 

2.  Take  the  column  of  Q  with  the  largest  number  of 
ones.   Assume  the  onqs  are  located  in  rows 
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3. 


4. 


i,j,  ...,k.   Form  the  submatrix  Y.  .      ,  by  re- 
taining  the  ith,  jth,  ...  kth  rows  and  columns  of 
Y  and  deleting  the  others.   Place  Y.  .     ,  in  a 
uniformly  tapered  form.   This  is  done  by  noting 
the  following 

(a)  Interchanging  row  and  column  j  with  row  and 
column  m  in  Y.  .     ,  corresponds  to  inter- 
changing  the  labeling  of  tree  branches  j  and 
m. 

(b)  Multiplying  all  elements  of  row  j  and  column 

j  by  (-1)  means  reversing  the  polarity  of  the 
jth  tree  branch ,    i.e.  the  jth  port.   (Note 
the  element  y. .  will  remain  positive). 
Repeat  step  2  for  the  columns  of  Q  with  successive- 
ly smaller  number  of  ones. 
Group  the  tree  branches  together. 


Example:   Realize  the  given  cut  set  matrix  Q  using  Kim's 
method. 


Q   = 


1 

3 

4 

e 


2 

5 

7 

8    i 

1 

3 

4 

6 

-1 

-1 

0 

! 

0    i 

1 

0 

0 

0 

1 

0 

0 

-1  1 

0 

1 

0 

0 

0 

0 

1 

1 

1 

0 

0 

1 

0 

1 

1 

1 

0 

0 

0 

0 

1 

QQ  =Y  = 


1 

3 

4 

6 

1 

3 

-1 

0 

-2 

3 

-1 

3 

-1 

1 

4 

0 

-1 

3 

1 

6 

-2 

1 

1 

4 
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Column  2  of  Q  has  the  maximum  number  of  non-zero  elements. 


■1,3,6 


13  6 

3   -1  -2 

-1    3  1 

-2    1  4 


This  matrix  is  not  uniform- 
ly tapered  since  it  violates 
condition  1,  i.e.  there  are 
negative  elements  in  Y. 


1,3,6' 


Multiplying  row  and  column  1  by  (-1)  yields 


\* 


1,3,6 


'I 

3 
6 


6 
,J 
1 

4 


This  matrix  is  still  not 
uniformly  tapered  since  it 
violates  conditions  2a  and 
2b. 
'Interchanging  rows  3  and  6  and  columns  3  and  6  (rule  a) 


1,6,3 


1 

6 

3 

1 

3 

2 

1 

6 

2 

4 

1 

3 

1 

1 

3 

which  is  uniformly  tapered, 


Proceeding  with  column  5  of  Q,  yields: 


1,6 


1 

6 

1 

3 

-2 

6 

-2 

4 

placing  in  uniformly 
tapered  form 


1 , 6 


Columns  7  and  8  yield 


1 
6 


2 

4 


Y 


6,4 


6 

4 

( 

1 

4 

1 

A 

and  Y 


3,4 


3 

3 

4 
-1 

placing  in 
uniformly 

3 

3 

4 

3 

3 

1 

4 

-1 

3 

tapered 
form 

4 

1 

3 
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The  linear  sub-trees  (Fig. 10)  are  then  combined  to  form  the 
tree  as  in  Fig.  9. 

It  can  be  seen   that  this  method  is  very  difficult   because 
the  process  of  obtaining  the  uniformly  tapered  form  is 
essentially  a  trial  and  error  procedure. 

3 . 3   Reducing  the  Cut  Set  Matrix  to  a  Set  of  Incidence 

Matrices 

Tutte  [ll]  and  Mayeda  [12]  derive  a  set  of  incidence 
matrices  from  a  given  cut  set  matrix  and  then  combine  their 
corresponding  subgraphs.   Auslander's  [13]  method  arrives 
at  similar  results  using  the  loop  concept  to  justify  the 

incidence  relations.   Tutte  bases  his  algorithm  on  the 

* 

theory  of  matroids  .   He  approaches  the  realization  problem 

in  much  the  same  way  as  Mayeda.   The  basic  ideas  underlying 

their  methods  are  as  follows: 

Connect  all  edges  of  a  fundamental  cut  set,  i  (Fig.  11a) , 

at  an  auxiliary  node,  n  (Fig. lib).   Then  split  this  node, 

thus  generating  two  separate  graphs  G,  and  G~ .   Removal  of 

the  cut  set,  i,  from  G,  and  G0  yields  two  subgraphs  G   and 

JL       &  a 

G,  .  (Fig.  lie) .   The  removal  of  the  cut  set,  i,  from  the 
graph  corresponds  to  omitting  row  i  from  Q  and  all  columns 

*  A  binary  matroid  is  the  class  of  elementary  chains  of  a 
binary  chain— group,  i.e.   the  edges  of  a  fundamental 
cut  set  form  an  elementary  chain,  all  fundamental  cut 
sets  with  respect  to  a  given  tree  form  a  binary  chain- 
group,  thus  constituting  a  class  or  binary  matroid.  This 
binary  matroid  can  be  represented  by  a  matrix,  which  is 
called  Q  in  this  thesis. 
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Figure  8 


Linear  tree  port  structure  of  a  graph. 


Figure  9 
Graph  of  example  in  Section  3.2. 
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Figure  10 
Subtrees  found  by  Kim's  procedure 
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of  Q  that  have  a  one  in  this  row.   This  operation  yields 
a  new  matrix  H.   If  H  can  be  partitioned  as 


H  = 


H.,^   0 


0    H, 


the  cut  set  i  does  not  represent  an  incidence  set,  that  is 
the  elements  of  i  do  not  converge  on  a  single  node.   Thus 
both  separate  parts  G   and  G,  contain  at  least  two  nodes, 
one  tree  branch,  and  one  link.   H,  and  H~  are  the  funda- 
mental cut  set  matrices  of  G   and  G,  ,  respectively. 
The  fundamental  cut  set  matrix  M, (i)  of  G,  is  found  by  re- 
moving from  Q  all  columns  and  rows  which  belong  to  H„ . 
Similarly  M~(i)  of  G~  is  found.   If  H  cannot  be  partitioned, 
one  of  the  separate  parts  has  to  consist  of  only  one  node, 
indicating  that  the  cut  set  i  is  an  incidence  set.   In 
this  case  M, (i)  =  Q  and  M~(i)  =  row  i,  omitting  the  zero 
entries.   Proceed  by  successively  removing  the  other  rows 
until  all  M-matrices  are  incidence  matrices.   If  Q  has  b 
rows,  (b+1)  of  these  "minimum  M-submatrices"  [12]  are 
obtained,  yielding  (b+1)  subgraphs.   These  subgraphs  are 
then  combined  by  joining  them  at  the  auxiliary  nodes. 
Example:   Realize  the  following  Q  matrix  using  Mayeda ' s 
procedure.   The  corresponding  graph  is  shown  in  Fig.  12. 


Q   = 


5 
6 

7 

8 


1 

2 

3 

4    i 
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6 

7 

8 

1 

1 

0 

0    ■ 

1 

0 

0 

0 

1 

I 

1 

1 

0 

1 

0 

0 

0 

1 

1 

1 

0 

0 

1 

0 

0 

1 

1 

0 

0 

0 

1 

0 
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Underlining  indicates  minimum  matrix. 


Remove  cut  set  i  =  5  from  Q 


H   = 


3 

4 

6 

7 

8 

1 

1 

1 

0 

1 

1 

1 

0 

1 

0 

1 

0 

0 

0 

1 

Partitioning  is  not  possible. 

12   5 
Letting  M1(5)  =  Q  and  M2(5)  =  5  |T~  1   1  ,  M2(5)  represents 

a  minimum  M-matrix. 

Proceeding  with  M, (5)  cut  set  j  -  6  is  removed. 


H   = 


5 

7 

8 

5 

1 

0 

0 

7 

0 

1 

0 

8 

0 

0 

1 

Omitting  rows  and  columns  7  and  8  in  Q  yields 


M,  (56)  = 


1 

2 

3 

4 

5 

6 

1 

1 

0 

0 

1 

0 

11110   1 


Omitting  row  and  column  5  in  Q  yields: 


M1(6) 


1 

2 

3 

4 

6 

7 

8 

1 

1 

1 

1 

1 

0 

0 

0 

1 

1 

1 

0 

1 

0 

Proceeding  with  M-,  (6)  cut  set  k  =  7  is  removed. 
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k  =  7 


Mn(67) 


1 

2 

3 

4 

6 

7 

1 

1 

1 

1 

1 

1 

0 

1 

1 

1 

0 

1 

M1(7) 


2 

3 

4 

7 

8 

7 

1 

1 

1 

1 

0 

8 

1 

1 

0 

0 

1 

Finally  cut  set  1  =  8  is  removed  from  M, (6) 


=  8 


M1;L(68)  =   M1(6) 


1^(8) 


2   3   8 
=   8  I  1   1   1 


The  realization  steps  are  shown  graphically  in  Fig.  13. 


3 .4   Successive  Removal  of  Tree  Tips 

The  methods  of  Guillemin  [14]  and  Iri  [15]  successive- 
ly remove  the  outermost  branches  of  a  tree  ( " tips" [14]) , 
once  these  have  been  identified.   Since  Iri's  method  is 
more  general  (Guillemin1 s  procedure  fails  in  the  case  of 
star  trees)  its  principles  shall  be  outlined.   The  main 
idea  is   that  only  those  edges  are  removed  which  do  not 
make  the  resulting  graph  separable.   This  is  achieved  by 
three  operations. 

Thinning  (Fig.  14)  which  means  to  delete  a  column  with  only 
one  non-zero  element  or  to  delete  all  but  one  of  a  set  of 
identical  columns.   This  removes  parallel  links. 

3      4      5  3 


-F 


110 
111 


becomes 
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Shrinking  (Fig.  15) ,  which  means  to  delete  rows  (1  and  2) , 
which  have  only  one  non-zero  element  in  the  same  column. 
If  there  are  identical  rows  (4,6,7),  one  interchanges  the 
role  of  a  link  (10)  and  a  tree  branch  (4) ,  corresponding 
to  one  of  these  rows.   The  result  is  that  all  but  one  of 
these  rows  have  only  one  non-zero  element  in  a  common 
column  (column  4) .   The  branches  corresponding  to  these 
rows  are  "perfectly  series"  [15],  i.e.  they  are  in  series 
in  any  2-isomorphic  graph  (cf.  sect.  3.5) 


1 
2 
3 
-Ffc  =  4 
5 
6 
7 


8 

9 

10 

1 

0 

0 

1 

0 

0 

0 

0 

1 

0 

1 

1 

1 

1 

1 

0 

1 

1 

0 

1 

1 

8 

9 

4 

1 

1 

0 

0 

2 

1 

0 

0 

3 

0 

1 

1 

10 

0 

1 

1 

5 

1 

1 

1 

6 

0 

0 

1 

7 

0 

0 

1 

8   9   4 


f 
2 

1 

0 

0 

3 

0 

1 

1 

0 

0 

1 

1 

5 

1 

1 

1 

6 

0 

0 

1 

The  left  hand  form  of  -F  corresponds  to  Fig.  15a,  the 
right  hand  form  to  Fig.  15c,  the  one  in  the  center  to  Fig. 
15b. 

Reduction  (Fig.  16)  which  means  to  remove  a  "perfect  tip" 
[15],  i.e.  a  tree  branch,  which  is  a  tip  in  any  2-isomorphic 
graph.   The  reduction  with  respect  to  row  1  is  obtained  as 
follows:   add  or  subtract  pairs  of  those  columns  (6  and  9) 
having  a  non-zero  element  in  row  1,  such  that  the  element 
of  the  resulting  column  has  a  zero  in  row  1.   Augment  the 
matrix  by  those  of  the  new  columns,  which  have  more  than 
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1   ' 


(a) 


- 


y 


(b) 


C   £*> 


(c) 


Fig.  11  -  Reduction  of  a  graph  by 
coalescing  and  removing 
a  cut  set. 


Fig. 12  - 


Graph  of 
example  in 
Section  3.3 
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H     '56) 


Mu  (e>7)  M>  (?) 
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Figure  13 
Realization  steps  of  Mayeda's  procedure 
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one  non-zero  element  and  which  are  not  identical  with  any 
column  already  present  in  the  matrix  (column  a) .  Finally 
remove  row  1.  The  reduction  procedure  attaches  to  a  graph 
new  links  (a)  corresponding  to  those  columns  by  which  the 
matrix  is  augmented.  (Each  of  these  links  forms  a  3  edge 
loop  together  with  the  two  links  (6  and  9)  connected  to  a 
perfect  tip  1) .  Removal  of  row  1  corresponds  to  removing 
this  perfect  tip.  (Fig.  16b) 


8 


-F 


1 
2 

3 
4 

5 


10  0  10 
10  110 
-1-1  0  0  1 
0  1-1-1  -1 
0-1100 


Example:       Realize   Q   =    [-F      u]    by   Iri's   method.       This   cut 
set  corresponds    to    the   graph   shown   in   Fig.    17a, 


-F      = 


1 

2 

3 

4 

-1 

1 

0 

0 

1 

-1 

-1 

-1 

0 

1 

1 

1 

0 

-1 

-1 

0 

No  thinning  or  shrinking  are  possible.   Thus  reduce  with 
respect  to  row  5 


1 

2 

3 

4 

a 

5 

-1 

1 

0 

0 

0 

6 

1 

-1 

-1 

-1 

0 

7 

0 

1 

1 

1 

1 

8 

0 

-1 

-1 

0 

-1 
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3      ¥ 


(a) 


(b) 


Figure  14 
Illustration  of  "Thinning" 
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10 


(a) 


(b) 


^  ,   »   ,  6  v,   3 


Figure  15 
Illustration  of  "Shrinking" 


\8 


Zf 
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:     \ 

<. 

»V7 

Figure  16 
Illustration  of  "Reduction" 
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6 


1 

2 

3 

4 

a 

1 

-1 

-1 

-1 

0 

0 

1 

1 

1 

1 

0 

-1 

-1 

0 

-1 

Thin  with  respect  to  column  1 


2   3   4 


7 

1 

1 

1 

1 

8 

-1 
2 

-1 

0 

-1 

7 

1 

8 

-1 

7   0   1111  then  augment  and  reduce  with 

respect  to  row  6. 

Thin  with  respect  to  column  4, 
then  remove  all  but  column  2. 

Shrink  with  respect  to  rows  7 
and  8. 

The  result  is  a  0  x  1  "void"  [14]  matrix,  representing  a 

single-edge  loop. 

The  graph  is  realized  starting  with  the  single-edge  loop. 

The  links  and  tree  branches  are  added  to  the  graph,  which 

have  been  removed  by  the  operations  performed  on  the  cut 

set  matrix  (Figs.  17b  to  17h) . 

This  procedure  rests  heavily  on  theorems  whose  proof  is  in 

some  cases  "of  considerable  length  and  tedious"  [15 ] . 

3 . 5   Combining  Loops  within  2-isomorphisms 

Two  concepts  have  to  be  introduced  before  the  methods 
of  Fu  [16]  and  Lofgren  [17]  can  be  discussed. 

1.   In  a  switching  network  each  branch  corresponds  to 
a  switch,  which  is  closed  (1)  or  open  (0) .   Con- 
necting two  nodes  by  parallel  paths  generates  a 
switching  function.   If  one  switch  within  a  path 

*   A  single-edge  loop  is  formed  by  coalescing  the  two  nodes 
of  an  edge. 
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is  open,  the  whole  path  is  open.   If  one  path  is 
"closed",  the  two  nodes  A  and  B  are  connected 
(Fig.  18) . 
2.   Two  graphs  G,  and  G~  are  2-isomorphic  if  G„  can  be 
found  from  G,  by  following  operations  (Fig.  19a) : 
Cut  the  graph  G,  at  two  nodes  (n  and  m)  into  dis- 
joint subgraphs.   Turn  one  of  them  around  and  join 
them  again,  yielding  G~.   Two  graphs  are  1-isomor- 
phic   if  they  are  cut  or  joined  at  one  node  (Fig. 
19b). 
Fu's  method  is  based  on  the  following  observation:   Short 
circuiting  a  tree  branch  j  (deleting  column  j  in  F) ,  yields 
a  new  graph  G . ,   which  has  the  same  loops  as  the  original 
graph  G,  except  for  those  loops,  which  contain  the  short 
circuited  branch  j.   In  Figs.  20a  and  20c  the  tree  branch  2 
is  short  circuited.   If  the  resulting  graph  G-  is  non-sep- 
arated (Fig.  20b) ,  the  loops  which  contained  branch  2  in 
the  original  graph  will  not  contain  branch  2  in  G?.   If  G~ 
is  separable  (Fig.  20d) ,  the  loops  which  contained  branch 
2  in  the  original  graph  will  be  single-edge  loops.   The 
underlying  idea  of  LOfgren's  method  is  that  tree  branches 
within  each  loops  have  to  form  an  uninterrupted  linear  path 
and  that  the  order  of  tree  branches  common  to  more  than  one 
loop  has  to  be  the  same  within  each  loop.   During  the  real- 
ization procedure  those  branches  common  to  two  loops  are 
rearranged  within  2-isomorphisms.   This  is  illustrated  in 
Fig.  21.   Two  loops  (I  and  II)  shall  be  combined.   Three 
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(a) 


O 

(b) 


V  em 


(c) 


(d) 


(e) 


(f) 


(9)  (h) 

Figure  17 
Realization  steps  in  Iri's  procedure. 


Figure  18 
Switching  network 
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branches  1  and  3  do  not  form  an  uninterrupted  path  in  loop 
I.   Therefore  a  2-isomorphic  form  of  loop  I  is  constructed 
(Fig.  21c) ,  in  which  branches  1  and  3  form  a  linear  path 
before  the  loops  are  joined  (Fig.  21d) . 
The  realization  steps  in  Fu ' s  method  are: 

(a)  Short  circuit  all  but  one  of  the  tree  branches 
and  draw  the  corresponding  graph. 

(b)  Successively  insert  the  other  tree  branches  in 
the  graph  following  the  conditions  demonstrated 
in  Fig.  20. 

Although  this  method  looks  simple  it  is  not  systematic  be- 
cause a  tree  branch  can  often  be  inserted  in  more  than  one 
place.   This  yields  alternate  subgraphs,  some  of  which  will 
not  lead  to  a  valid  graph.   These  "dead  graphs"  [16]  have 
to  be  carried  along  to  ensure  that  a  valid  graph  will  result, 
The  realization  steps  of  Lofgren's  method  are: 

(a)  Choose  the  loop  which  contains  the  smallest  number 
of  tree  branches  and  draw  the  corresponding  graph. 

(b)  Expand  this  graph  by  those  loops  which  share  the 
smallest  number  of  tree  branches  with  the  already 
existing  graph.   If  these  common  branches  do  not 
form  an  uninterrupted  linear  path,  rearrange  them 
within  2 -isomorphisms. 
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(a)  (b) 

Figure  19 
2-isomorphic  (a)  and  1-isomorphic  graphs  (b) . 


/      2 
* ^ 


(a) 


(c) 


(b) 


i 

(d) 


Figure  20 


Non  separated  (b)  and  separable  graph  (d)  resulting  from 
removal  of  branch  from  original  graphs  (a  and  c) . 


3 

-• » 


vJ    Z 


_  H-  - *_ 

(a)  (b)  (c)  (d) 

Figure  21 
Combining  loops  (a,b)  within  2 -isomorphisms  (c)  to  a  graph  (d) 
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1 

2 

3 

4 

5 

6 

7 

8 

1 

0 

0 

0 

1 

1 
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0 
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0 

0 

1 

1 

1 

1 

0 

0 

1 

0 

0 

1 

1 

1 
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1 

0 

1 

1 

0 

Example:   Realize  B  using  Fu ' s  and  Lof gren ' s  methods.   The 
given  loop  matrix  corresponds  to  the  graph  in 
Fig.  12. 


1 

B   =     2 
3 

4 


The  realization  steps  are  shown  in  Fig.  22. 

3 .6   Generating  the  Tree  from  Path  Sets 

This  method,  due  to  Gould  [18],  treats  the  tree 
branches  of  each  loop  as  a  path  set.   The  tree  branches 
common  to  different  loops  form  additional  path  sets.   Or- 
dering maximal  path  sets,  i.e.  sets,  which  are  not  contain- 
ed in  any  other  path  set,  subtrees  are  obtained.   This 
ordering  procedure  is  done  by  using  those  path  sets  which 
are  contained  in  the  respective  maximal  sets.   This  proce- 
dure is  extremely  laborious  if  the  tree  is  not  a  linear 
path. 
Example:   Demonstrate  Gould's  procedure.   The  loop  set 

matrix  corresponds  to  Fig.  23a.   The  graphical 
representation  of  the  maximal  path  sets  and  the 
forming  of  the  tree  is  given  in  Fig.  23b. 


49 


B       = 


1 

2 
4 
8 


12483567 

10000110* 

01001100* 

00100011 

00011011* 


set   1   and   2:     (3(5)6)  *denote 

set  4    and   8:     (3(67))  maximal 

path   sets 

These  sets  are  combined  using  "overlap  rules"  [18],  which 
merely  state  that  one  should  combine  the  subtrees  by  over- 
lapping corresponding  branches  to  obtain  the  whole  tree. 

3.7   Summary 

The  basic  principles  of  existing  procedures  have  been 
classified  and  outlined  with  the  most  efficient  method  of 
each  classification  used  to  solve  a  simple  realization  pro- 
blem.  Each  of  these  methods,  however,  suffers  from  at 
least  one  of  the  following  deficiencies  in  its  application 
as  an  engineering  tool. 

1.  The  method  fails  in  some  cases  (Guillemin) . 

2.  The  underlying  mathematical  principles  are  too 
abstract  for  practical  applications  and  thus  do 
not  give  much  insight  in  the  realization  proce- 
dure (Tutte,  Iri) . 

3.  Except  for  Iri's  method,  none  of  these  procedures 
has  been  put  into  the  form  of  an  algorithm,  which 
can  be  directly  used  on  a  digital  computer. 
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Figure  22 
Realization  steps  in  Fu ' s  (a)  and  LOfgren's  method  (b) 


,   3   l £ . 


J  6       7 

•       •  •       » 

S 
»       \i  »■  -  ■  « 


(a)  (b) 

Figure  23 
Demonstration  of  Gould's  procedure 
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4.   The  method  requires  too  much  work  particularly  as 
the  matrices  become  more  complicated.   This  is 
true  for  all  procedures  reviewed. 

The  last  statement  is  especially  true  for  those  methods 
which  reveal  alternatives  while  the  realization  progresses, 
since  only  some  of  these  alternatives  will  yield  a  valid 
graph.   Thus  one  has  to  carry  along  "alternative"  subgraphs, 
which  finally  may  end  up  as  being  a  "dead  graph"  [16], 
Moreover,  rewriting  matrices  and  redrawing  graphs  is  always 
a  source  of  error. 
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REALIZATION  OF  THE  LOOP  MATRIX  B 
USING  THE  COMBINATORIAL  APPROACH 

4.1  Introduction 

The  method  presented  in  this  thesis  takes  a  combina- 
torial approach,  i.e.  joining  I    loops  together  such  that 
common  branches  appear  in  uninterrupted  chains.   It  thus 
belongs  to  the  classification  of  "Generating  the  Tree  from 
Path-Sets".   This  method  identifies  trunk  branches,  main 
branches,  and  limbs,  and  introduces  the  concept  of  "Unique 
connections".   This  yields  a  procedure  which  is  easy  to 
apply  and  generally  more  efficient  than  those  reviewed 
above.   The  required  theorems  to  develop  the  procedure  are 
proven  rigorously  using  the  method  of  contradiction.   Also 
the  method  has  been  programmed  for  automatic  computation  on 
a  digital  computer  as  discussed  in  Chapter  5. 

4. 2  The  Concept  of  Trunk  Branches,  Main  Branches,  and  Limbs 
The  following  information  about  a  graph  associated  with 

a  given  loop  matrix  can  be  obtained  immediately  from  the 
B-matrix: 

(a)  The  number  of  loops  in  the  graph,  and  the  tree 
branches  plus  the  defining  link  for  each  loop. 

(b)  Those  tree  branches  which  are  common  to  more  than 
one  loop. 

To  extract  more  information  five  theorems  are  stated. 

Theorem  4 ;   In  each  loop  the  k  tree  branches  form  a  linear 
path  whose  end  nodes  are  connected  by  a  link. 
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Proof:       Assuming  there  existed  a  loop  whose  tree 

branches  do  not  form  a  linear  path.   At  least 
one  tree  branch  would  share  a  node  with  two 
branches  in  a  linear  path.   This  would  contra- 
dict the  definition  of  a  loop  requiring  that 
exactly  two  edges  have  to  be  incident  with 
each  node. 

Property  1;   If  k  tree  branches  belong  to  only  one  loop, 

they  may  appear  within  the  linear  path  in  any 
order. 

This  property  holds  because  permuting  any  two  tree  branches 
does  not  violate  the  definition  of  a  loop. 

Theorem  5 ;    Two  or  more  tree  branches  common  to  more  than 

one  loop  have  to  appear  in  the  same  sequence  in 
all  loops. 

Proof:       Assuming  in  Fig.  24  tree  branches  (b,c,d)  form 
a  linear  path  b-c-d  in  loop  i  and  b-d-c  in 
loop  j.   This  is  possible  only  if  c  and  d  are 
identical,  i.e.  tree  branches  b  and  c  form  a 
linear  path  in  both  loops. 

As  a  result  of  Property  1  and  Theorem  5  only  those  tree 
branches  common  to  more  than  one  loop  will  put  restrictions 
on  the  ordering  of  tree  branches  within  loops.   Those 
branches  that  appear  in  only  one  loop  (their  corresponding 
column  in  the  B-matrix  has  only  one  non-zero  element) ,  can 
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be  placed  in  series  with  the  link.   We  will  call  them  limb 

* 

branches,   or  just  limbs   (branch  9  in  Fig.  25) .   The  loop 

matrix,  B  of  the  graph  in  Fig.  25  is 


B  = 


1 
2 
3 


1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

0 

0 

1 

1 

1 

1 

0 

0 

0 

1 

0 

1 

1 

0 

0 

1 

1 

0 

0 

1 

0 

0 

1 

1 

1 

0 

links 


tree  branches 


That  loop  which  contains  the  maximum  number  of  tree  branches, 
which  are  not  limbs,  shall  be  called  main  loop  (loop  1  in 
Fig.  25).   These  tree  branches  (branches  4,5,6,7)  form  the 
trunk*  of  the  graph,  and  are  called  trunk  branches'* 

Theorem  6 :   Trunk  branches  form  a  linear  path  not  interrupt- 
ed by  main  branches  or  limbs. 

Proof:       If  there  were  a  limb  within  the  trunk,  it  could 
always  be  moved  to  one  of  the  two  end  nodes 
since  it  appears  only  in  one  loop.   If  there 
were  a  main  branch  present  then,  by  definition, 
it  should  have  been  included  in  the  trunk. 

4 . 3   Ordering  Trunk  Branches  and  the  Concept  of  Unique 

Connections 

By  Theorem  5  each  group  of  trunk  branches,  belonging 
to  one  loop,  also  appears  in  the  trunk,  not  interrupted  by 


*   These  definitions  are  original.   The  term  "limb"  appears 
sometimes  in  the  literature  but  with  a  different  meaning, 
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any  other  branch.   These  groups  shall  be  called  chains  of 
trunk  branches  or  just  chains .   A  main  branch,  together 
with  some  of  the  trunk  branches,  form  linear  paths  in  at 
least  two  different  loops.   These  groups  of  trunk  branches 
also  form  chains  and  must  share  a  common  node  (cf.  Theorem 
6)  to  which  the  main  branch  is  connected  (Fig.  25,  node  a) . 
Thus  the  following  property  is  obtained: 

Property  2 :   The  chain  and  the  common  node  requirements  are 
the  only  restrictions  on  the  ordering  of  trunk 
branches . 
Once  the  trunk  is  established  by  combining  the  chains, 
starting  with  the  smallest  ones,  the  common  node  require- 
ments can  again  be  used  to  attach  the  main  branches  to  the 
trunk. 

To  ensure  an  efficient  connection  procedure,  those  main 
branches  are  attached  first  whose  position  relative  to  the 
trunk  is  completely  specified,  i.e.  no  other  branch  could 
take  this  location  without  causing  an  unrealizable  graph. 

Theorem  7a:   A  main  branch  m  is  uniquely  connected  to  the 
common  node  between  adjacent  trunk  branches  a 
and  b,  if  either  a  or  b  is  in  each  loop  in 
which  m  is  present,  provided  m  appears  at 
least  once  with  a  or  b,  respectively. 

Proof:        Assume  m  forms  a  linear  path  with  trunk  chains 
Cl  and  C2  in  loops  Ll  and  L2 .   If  CI  and  C2 
would  not  share  a  common  node,  they  must  be 
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separated  by  at  least  one  trunk  branch. 
Since  m  forms  a  linear  path  with  Cl  and  C2, 
respectively,  m  has  to  be  connected  to  an  end 
node  of  Cl  and  C2,  thus  forming  a  loop  (Fig. 
26a) ,  which  is  not  possible,  since  trunk  and 
main  branches  are  part  of  a  tree. 
Corollary  7;  If  more  than  one  main  branch  is  uniquely  con- 
nected to  the  same  node  and  there  exists  a 
loop  in  which  only  one  of  these  branches  ap- 
pears, this  branch  should  be  connected  first. 
If  there  is  no  loop  containing  all  of  these 
main  branches,  indicating  disjoint  sets  of 
main  branches,  one  branch  of  each  set  will  be 
connected  to  the  common  node. 

Theorem  7b :   A  main  branch  m  is  uniquely  connected  to  the 
end  node  of  a  chain,  if  the  trunk  branch  t, 
incident  with  the  end  node  appears  in  each 
loop  in  which  m  is  present,  provided  the  trunk 
branch  t~  adjacent  to  t..  appears  in  one  but 
not  all  of  these  loops. 

Proof:       Assume  m  were  connected  between  t,  and  t_, 

then  m,  t.  and  t„  could  not  appear  in  the  same 
loop.   Assume  m  were  connected  to  the  node  of 
t~,  at  which  t-  is  not  incident,  then  there 
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exists  a  loop,  in  which  m  has  to  be  parallel 
to  t2  (Fig.  26b) .  This  is  impossible,  since 
m,  t1  ,  and  t„  are  part  of  a  tree. 

Theorem  7b  will  also  be  used  to  attach  main  branches  to 
those  already  connected  to  the  trunk. 

Once  all  the  unique  connections  have  been  carried  out, 
main  branches  are  added  to  the  existing  tree  structure  by 
using  Theorem  4,  i.e.  they  have  to  form  linear  paths  with 
the  branches  already  connected. 

Should  the  graph,  corresponding  to  the  loop  matrix  B, 
be  separable,  this  will  be  indicated  by  the  following  pro- 
perty. 
Property  3 ;   If  there  exist  main  branches  or  limbs  which 

do  not  appear  in  any  loop  together  with  trunk 
branches,  they  belong  to  separate  parts  of 
the  graph,  none  of  which  contains  trunk  branch. 

In  case  there  should  be  no  main  branches  and/or  limbs,  the 
realization  procedure  will  not  be  affected  (see  examples  in 
sections  4.5.1  and  4.5.2). 
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Figure  24 
Illustration  of  Theorem  5 
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Figure  25 

Trunk  branches  (4,5,6,7),  main  branches  (8), 
and  limbs  (9) . 
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Figure  26 
Illustration  of  Theorems  7a  and  7b 
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4 .4   An  Illustrative  Example 

The  procedure  shall  be  demonstrated  by  working  out  a 
realization  in  full  detail. 
Given  the  matrix   B  =  [iJ'F],  where: 


f 


F  = 


Links 
(Loops) 


<— 
1 

2 

3 

4 

Tree 
5   6 

Branch 
7   8 

es  - 
9 

10 

11 

12 

-> 

13 

14 

0 

1 

0 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

15 

0 

1 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

16 

0 

0 

0 

0 

1 

0 

0 

1 

0 

1 

0 

0 

0 

17 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

18 

1 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0 

19 

0 

1 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0 

20 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

21 

1 

1 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0 

1 

22 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

23 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

1 

24 

0 

1 

1 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

(1) 


Since  the  procedure  does  not  involve  the  assignment  of  di- 
rections to  the  branches,  the  F  matrix,  as  given,  is  non- 
oriented.   Should  an  oriented  graph  be  required  the  appro- 
priate arrows  may  be  placed  on  the  branches  afterwards  to 
correspond  to  the  +  or  -  directions  associated  with  each 
unity  element  in  the  given  circuit  matrix. 

4.4.1   Partitioning  the  F  Matrix 

(a)   From  F  form  a  submatrix,  F<,  of  those  columns 

which  have  only  one  non-zero  element.   This  se- 
parates the  limbs.  The  ordering  of  columns  with- 
in F.    is  not  important.   The  submatrix  containing 
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the  remaining  columns  is  designated  by  F 


tm* 


Thus 


=  Cf 


tm 


F  ] 


For  the  given  example  the  limbs  correspond  to  tree  branches 
10,  12,  and  5  and  F  is  partitioned  as  follows: 


F  = 


1 

2 

3 

4 

6 

7 

8 

9 

13 

11 

10 

12  5 

14 

0 

1 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0  0 

15 

0 

1 

1 

0 

0 

1 

0 

0 

0 

1 

0 

0  0 

16 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0  1 

17 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0  0 

18 

1 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0  0 

19 

0 

1 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0  0 

20 

0 

1 

1 

0 

0 

0 

0 

0 

0 

1 

0 

0  0 

21 

1 

1 

0 

1 

0 

0 

0 

1 

1 

0 

0 

0  0 

22 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

1  0 

23 

0 

0 

0 

1 

0 

0 

0 

0 

1 

0 

0 

0  0 

24 

0 

1 

1 

0 

0 

0 

0 

1 

0 

1 

0 

0  0 

tm 


->  «e 


•pr 


(b)   Find  the  row  of  F.   which  contains  the  maximum 

x  tm 

number  of  non-zero  elements.   Separate  these 
columns  and  form  a  submatrix  designated  by  F. . 
The  branches  corresponding  to  these  columns  are 
the  trunk  branches.   The  ordering  of  columns  F, 


is  not  important.   The  remaining  columns  of  F 


tm 


form  the  submatrix  F  .   The  branches  correspond- 

m  c 

ing    to    these   columns  within   F      is   not   important, 
Partition   F   as    follows 


F      =       [P.'F     iF,] 
1 1    m ,     a 
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For  the  example ,  the  row  corresponding  to  loop  21  is  seen  to 
contain  five  branches  of  F.  ,  namely  1,2,4,9,  and  13; 
which  are  now  identified  as  trunk  branches.   In  those 
cases  where  more  than  one  row  contains  the  maximum  number 
of  non-zero  elements,  alternate  selections  of  trunk 
branches  are  possible.  For  the  example.  F,  is  now  partitioned 
as  follows: 


F   = 


1 

2 

4 

9 

13 

3 

6 

7 

8 

11 

10 

12 

5 

14 

0 

1 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

15 

0 

1 

0 

0 

0 

1 

0 

1 

0 

1 

0 

0 

0 

16 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

1 

0 

1 

17 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

0 

0 

18 

1 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

19 

0 

1 

0 

0 

0 

0 

1 

0 

1 

0 

0 

0 

0 

20 

0 

1 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

21 

1 

1 

1 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

22 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

1 

0 

23 

0 

0 

1 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

24 

0 

1 

0 

1 

0 

1 

0 

0 

0 

1 

0 

0 

0 

->  -<r- 


m 


-*>  «^- 


4.4.2   Establishing  the  Trunk 

From  the  F  matrix  list  all  of  the  trunk  branches  which 
are  contained  in  each  row,  for  those  rows  that  contain  more 
than  one  trunk,  branch.   Thus  for  the  given  example 


Loop 

14  and  24 

21 

23 


Trunk  Branches 

2-9 
1-2-4-9-13 

4-13 


Table  1 
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From  F.  for  each  main  branch,  list  those  trunk 
m 

branches  which  appear  in  the  same  row.   When  the  main 
branch  is  contained  in  more  than  one  row  write  these  trunk 
branches  on  the  same  line  and  separate  these  numbers  from 
the  others  by  a  comma.   For  the  given  example  this  listing 
becomes : 


Main  Branch 

Trunk  Branch 

Chains 

3 

2, 

2-9 

6 

2-9, 

1,  2 

7 

2, 

9 

8 

1, 

2 

11 

2, 

2-9 

Table  2 

Main  branch  6  is  connected  to  trunk  branches  2  and  9 
in  loop  14.   Branch  6  is  also  connected  to  trunk  branch  1 
in  loop  18,  and  to  trunk  branch  2  in  loop  19.   Commas  are 
used  to  separate  combinations  or  chains  of  trunk  branches 
which  have  the  same  main  branch  in  common.   These  chains  of 
trunk  branches  must  have  at  least  one  node  in  common  (Theo- 
rems 6  and  7) .   The  order  in  which  branches  are  listed 
within  a  chain  by  itself  is  not  important,  since  we  look  at 
only  one  loop  at  a  time  (Property  1) . 

Using  the  second  table,  starting  with  rows  involving 
the  smallest  number  of  chains  of  trunk  branches,  arrange 
the  trunk  branches  in  a  sequence  which  satisfies  all  of  the 
common  node  requirements.   Thus,  for  the  given  example, 
trunk  branch  chains  1,2  and  2,9,  respectively,  must  have  a 
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common  node.  This  establishes  the  trunk  branch  chain  1-2-9. 
The  tabulated  requirements  that  branch  2  and  chain  2-9  have 
common  nodes,  and  that  branches  1  and  2  and  chain  2-9  have 
common  nodes  is  also  satisfied  automatically.  Should  any  of 
the  common  node  requirements  be  contradictory  the  graph 
would  not  be  realizable. 

Inspect  this  sequence  and  add  additional  trunk  branches 
so  that  all  of  the  trunk  branches  in  Table  1  are  included  in 
the  trunk  chain.  This  establishes  the  trunk  of  the  tree.  For 
the  example,  trunk  branches  4  and  13  must  be  added.  Since  4 
and  13  must  be  in  a  chain  it  follows  that  the  trunk  may  be 
established  in  four  alternate  ways  as  shown  in  Fig. 27,  the 
only  restrictions  being  that  chains  1-2-9  and  4-13,  respect- 
ively, be  kept  intact. 

4.4.3   Making  Unique  Connections 

The  information  contained  in  Table  2  is  now  used  to 
carry  out  unique  connections,  applying  Theorems  7a,  7b,  and 
Corollary  7. 

Main  branch  7  is  uniquely  connected  (Theorem  7a)  to  the 
node  shared  by  trunk  branches  2  and  3. 

Main  branch  8  is  uniquely  connected  (Theorem  7a)  to  the 
node  between  1  and  2,  but  so  are  main  branch  6  (Theorems  7a 
and  7b) ,  main  branch  3  (Theorem  7b) ,  and  main  branch  11 
(Theorem  7b) .   Applying  Corollary  7  the  partitioned  F-matrix 
is  examined.   There  is  a  loop  (14)  containing  only  main 
branch  6  in  addition  to  trunk  branches,  whereas  loops  18 
and  19  contain  main  branches  6  and  8.   Thus  main  branch  6 
has  to  be  connected  first  and  is  followed  by  8. 

64 


Main  branches  3  and  11  do  not  appear  together  with  6 
or  8  in  any  loop,  thus  they  will  be  part  of  a  different 
linear  path  (Corollary  7) .   Since  there  is  no  loop,  which 
besides  trunk  branches  contains  only  3  or  11,  either  one  of 
them  may  be  connected  first.   Choosing  main  branch  3  the 
unique  connections  are  shown  in  Fig.  28. 

4.4.4   Extending  Main  Branches 

Continue  the  tabulation  by  listing  unconnected  main 
branches  and  all  chains  of  trunk  branches  and  connected 
main  branches  which  appear  in  the  same  row,  using  commas 
to  separate  chains  for  different  rows.   Inspect  this  part 
of  the  table  to  connect  additional  main  branches  to  the 
graph  as  in  section  4.4.3.   After  these  branches  are  con- 
nected extend  the  table  further,  continuing  until  all  main 
branches  are  connected  to  the  tree. 

For  the  given  example  the  remaining  unconnected  main 
branch  is  11. 
The  extended  tabulation  becomes; 

Main  Branch       Trunk  Branch  and  Connected  Branch 

Chains 

11  2-3-7,  2-3 

By  Theorem  7b  main  branch  11  is  uniquely  connected  to  main 
branch  3  at  the  node  where  trunk  branch  2  is  not  incident. 
This  completes  the  connection  of  main  branches. 
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4.4.5  Connecting  the  Limbs 

Add  the  limbs  to  the  graph  so  that  they  are  grouped 
with  the  trunk  and  connected  main  branches  as  indicated  by 
the  rows  in  the  F  matrix. 

4.4.6  Connecting  the  Links 

The  position  of  the  links  relative  to  the  tree  are 
uniquely  determined  by  the  B  matrix. 

The  application  of  the  last  two  steps  to  the  given  ex- 
ample yields  the  final  graph  as  drawn  in  Fig.  29. 

4 . 5   Further  Examples 

4.5.1   Graph  used  in  Review  of  Existing  Realization  Proce- 
dure (No  Limbs  Present) 


F  = 


3 

5 

6 

7 

1 

0 

1 

1 

0 

2 

1 

1 

0 

0 

4 

0 

0 

1 

1 

8 

1 

0 

1 

1 

The  trunk  branches  can  immediately  be  identified  as;  3,6, 
7 ,  thus  5  is  a  main  branch. 

Table  1  Table  2 

6-7  5:  6,3 

3-6-7 

Trunk t      1    -   6-3 

Since  main  branch  5  is  uniquely  connected  to  the  node  be- 
tween trunk  branches  6  and  3,  the  resulting  graph  is  that 
of  Fig.  30. 
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Figure  27 
Possible  trunks  of  illustrating  example. 
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Figure  28 
Unique  connections  of  illustrating  example. 
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Figure  29 
Final  graph  of  illustrating  example. 
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4.5.2   Graph  with  Separate  Parts  (No  Main  Branches  Present) 


8 
9 


F   = 


10 
11 


1 

2 

3 

4 

5 

6 

7 

1 

1 

1 

0 

0 

0 

0 

0 

1 

1 

1 

0 

0 

0 

1 

1 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

1 

m 


Table  1 
1-2-3 
2-3 
1-2 

Trunk:    1-2-3 

Possible  connections  for  limbs  4  and  5  are  shown  in  Fig. 31a, 
The  separate  part  consists  of  limbs  6  and  7  and  link  11. 
The  final  graph  is  shown  in  Fig.  31b. 

4.5.3   Loop  Matrix  which  is  not  realizable 

This  example  of  an  unrealizable  loop  matrix  has  been 
investigated  by  many  authors  (for  references  of  [3]). 


F  = 


5 

6 

7 

1 

1 

1 

0 

2 

I 

0 

1 

3 

0 

1 

] 

4 

1 

1 

1 

Trunk  Branches  5,6,7 
Table  1 
5-6 
5-7 
6-7 
5-6-7 


The  only  way  these  chains  can  be  combined  is  shown  in  Fig, 
32.   This  violates  Theorem  6  which  states  that  trunk 
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(a) 

Figure  30 
Graph  of  Section  4.5.1. 
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Figure  31 
Graph  of  Section  4.5.2. 


(c) 

Figure  3  2 
Trunk  of  Section  4.5.3 
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branches  have  to  form  a  linear  path.  Moreover  Theorem  4 
is  violated.  Thus  this  loop  matrix  is  not  realizable  as 
a  graph. 

4.5.4   Graph  Consisting  of  Trunk  and  Main  Branches,  Limbs, 
and  Separate  Parts 
Given: 


F  = 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

13 

14 

15 

16 

17 

0 

1 

1 

0 

0 

1 

0 

0 

0 

1 

1 

0 

0 

0 

0 

1 

18 

1 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

0 

0 

0 

19 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

1 

20 

1 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

1 

21 

1 

0 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

1 

0 

0 

0 

22 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

23 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

24 

0 

0 

1 

0 

0 

1 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

25 

0 

0 

0 

0 

1 

1 

0 

0 

0 

0 

0 

1 

0 

1 

0 

0 

26 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

1 

0 

After  partitioning 


F  = 


1  9 

13 

16 

3  5 

6  7 

8 

12 

14  2 

4 

10 

11 

15 

17 

0  0 

0 

1 

1  0 

1  0 

0 

0 

0  1 

0 

1 

1 

0 

18 

1  0 

0 

0 

0  0 

0  1 

1 

0 

0  0 

0 

0 

0 

0 

19 

0  0 

0 

1 

0  0 

0  0 

1 

0 

0  0 

0 

0 

0 

0 

20 

1  1 

1 

1 

0  0 

0  0 

0 

0 

0  0 

0 

0 

0 

0 

21 

1  0 

1 

0 

0  0 

0  1 

1 

0 

0  0 

0 

0 

0 

0 

22 

0  1 

0 

0 

1  0 

0  0 

0 

0 

0  0 

0 

0 

0 

0 

23 

0  0 

0 

0 

0  1 

1  0 

0 

1 

0  0 

0 

0 

0 

0 

24 

0  0 

0 

0 

1  0 

1  0 

0 

1 

0  0 

0 

0 

0 

0 

25 

0  0 

0 

0 

0  1 

1  0 

0 

1 

1  0 

0 

0 

0 

0 

26 

0  0 

0 

0 

0  0 

0  0 

0 

0 

0  0 

1 

0 

0 

1 

m 
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Table  1: 

1-9-13-16 

1-13 


Table  2  : 


3 

16, 

9 

- 

_ 

- 

5 

- 

- 

6 

6-12 

6 

16 

3-16,    3 

- 

- 

7 

1. 

13-1 

1-8,    1-13-8 

- 

- 

8 

1, 

16,    13-1 

- 

- 

- 

12 

- 

3 

6-3, 

6 

- 

The  realization  steps  and  the  final  graph  are  depicted  in 
Fig.  33. 
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Figure  33 
Realization  steps  and  final  graph  of  Section  4.5.4 
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PROGRAMMING  THE  REALIZATION  PROCEDURE 
ON  A  DIGITAL  COMPUTER 

5. 1   Introduction 

In  addition  to  the  merits  listed  in  section  3.7  there 
is  another  advantage  of  the  procedure  presented  in  this 
thesis:   All  realization  steps  are  "final"  and  a  connection 
will  never  be  nullified.   Thus,  when  this  procedure  reveals 
an  inconsistancy,  the  loop  matrix  is  not  realizable. 


5. 2   Input/Output  Considerations 

The  inputs  for  a  computer  program  have  to  be  a  mini- 
mum.  For  this  program  this  is  achieved  by  reading  in  an 

augmented  F  matrix,  F^ .   F  is  augmented  by  one  row  at  the 

— — - - — ^  a 

top  whose  elements  denote  the  tree  branches  corresponding 
to  the  columns  of  F,  and  by  one  column  at  the  left  side, 
indicating  the  links  or  loops  corresponding  to  the  rows  of 
F: 

0_  (_tr.ee.  bra^n_hes_ 

t 

i 

F         n 
a 

Thus  the  input  for  the  realization  of  the  B  matrix  associa- 
ted with  the  graph  in  Fig.  34  is: 
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The  output  or  print  out  imposes  more  problems.   There  are 
essentially  three  possibilities: 

1.  Graphical  output,  i.e.  have  the  computer  draw  the 
graph. 

2.  Numerical  output,  i.e.  using  a  number  code  or  the 
incidence  matrix. 

3.  Language  print  out,  i.e.  giving  instructions  how 
to  draw  the  graph. 

Possibility  1   requires  the  use  of  an  additional  difficult 
output  program.   Possibility  2   requires  the  use  of  an  ef- 
ficient coding  system  or  the  generation  of  the  incidence 
matrix. 

In  both  cases  information  about  unique  connections  is 
lost.   To  use  a  language  print  out,  however,  is  quite  logi- 
cal for  this  procedure,  because  the  graph  can  be  drawn  in 
steps  once:   the  possible  trunks  are  listed  in  their  proper 
sequence;   the  nodes  to  which  main  branches  are  uniquely 
connected  are  given;   similarly  the  nodes  to  which  the  re- 
maining main  branches  and  limbs  are  connected  are  listed; 
finally,  the  two  nodes  to  which  each  link  is  connected  are 
stated. 

5. 3   Partitioning  the  F  Matrix  and  Establishing  the  Trunk 

The  element  of  each  column  of  the  submatrix  F  within  F^ 

a 

are  summed.   If  the  sum  equals  one,  the  column  corresponds 

to  a  limb,  and  the  column  (9)  of  F   is  stored  in  an  array, 

a 

LIMB.   These  columns  form  F-.   Then  the  elements  of  each 
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row  in  F  are  summed,  except  those  elements  which  appear  in 

a  column  identified  as  a  limb.   The  row  (13)  having  the 

largest  sum  is  identified  as  the  main  loop.   Those  columns 

(2,3,4,5)  having  a  non-zero  element  in  this  row  correspond 

to  trunk  branches  and  are  stored  in  the  array,  TRUNBR. 

These  columns  form  F  .   Those  columns  (6,7,8)  having  a  zero 

in  this  row  correspond  to  main  branches  and  are  stored  in 

the  array  MAINBR.   These  columns  form  F  .   If  there  are 

m 

more  than  one  row  with  the  same  maximum  sum,  the  uppermost 
row  of  them  is  identified  as  the  main  loop;  alternate  se- 
lections will  yield  the  same  graph.   This  completes  the 
partitioning  step. 

From  F   all  chains  of  trunk  branches  are  stored  in  the 
array  KCHAIN,  employing  row  1  of  F  .   Thus  5678,67,58  are 

a. 

stored.   From  F   and  F^  the  common  node  requirements  for 
m      t  ^ 

each  main  branch  are  listed  in  KCHAIN,  avoiding  duplica- 
tions.  Thus  57  and  67  are  stored.   The  next  step  is  to 
combine  these  chains.   This  combination  procedure  is  some- 
what intricate  and  shall  be  explained  in  more  detail. 
Starting  with  the  smallest  chains,  stored  in  the  array 
KCHAIN,  larger  chains  are  formed  by  combinations,  stored 
in  the  array  JTRCH.   To  ensure  that  the  elements  are  cor- 
rectly ordered,  the  following  technique  is  used:   If  there 
is  no  restriction  on  the  order  in  which  the  elements  of 
the  resulting  JTRCH  have  to  appear,  variations  are  generat- 
ed such  that  each  element  appears  at  least  once  as  the 
first  or  last  element.   If  there  is  a  subsequent  KCHAIN, 
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which  has  at  least  one  element  in  common  with  a  JTRCH, 
then  there  exists  one  variation  of  this  JTRCH  which  has 
this  element  in  an  end  position.   This  variation  has  the 
correct  order,  i.e.  it  appears  in  this  order  also  in  the 
final  trunk  or  trunks.   Since  the  KCHAIN  reflects  only  a 
common  node  requirement,  those  elements  of  KCHAIN  which 
are  not  present  in  this  JTRCH  are  attached  to  the  latter. 
The  elements  by  which  the  JTRCH  is  augmented  are  referred 
to  as  "non-common"  elements. 
Ex, :    KCHAIN    38 

JTRCH     1234,  2413 
new     JTRCH     24138 

Once  a  JTRCH  has  been  augmented,  all  variations  of  it  are 
deleted,  i.e.  1234  and  2413.   The  underlying  idea  of  this 
technique  is  that  trunk  branches  form  an  uninterrupted  path 
(Theorem  6,  sect.  4.2).   This  combination  procedure  is  es- 
sentially a  problem  of  pattern  recognition  which  is  straight 
forward  for  a  person  to  perform.   However,  the  computer  has 
to  be  programmed  to  distinguish  between  7  combination  pat- 
terns. 
If  a  KCHAIN  (K) 

1.  cannot  be  combined  with  any  of  the  JTRCH,  then  it 
forms  a  new  JTRCH  (Jl) .   Variations  (J2)  of  this 
new  JTRCH  are  generated  as  described  above. 

K:  1234;   Jl:  1234;   J2  2143 

2.  indicates  that  some  of  the  JTRCH  (J)  violate  a 
common  node  requirement,  then  delete  these  JTRCH's. 
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K:  234,   J:  2143 

3.  is  completely  contained  in  a  JTRCH  (Jl) ,  then  no 
combinations  with  K  are  carried  out,  but  those  of 
the  JTRCH' s  are  deleted  which  violate  a  common 
node  requirement  (J2) . 

K:  578,   Jl:  48572,   J2 :  84572 

4.  contains  the  first  (or  last)  element  of  a  JTRCH  (J), 
then  a  new  JTRCH  (JN)  is  formed,  augmenting  J  by 
the  "non-common"  elements  of  K.   All  variations  of 
J  are  deleted. 

K:  7453,   J;  1234,   JN:  123475 

The  variations  of  JN  are  123457,  thus  each  of  the 
new,  unordered  elements  is  once  in  an  outer  posi- 
tion. 

5.  contains  the  first  (or  the  last)  element  of  two 
different  JTRCH' s  (Jl  and  J2) ,  then  Jl  and  J2  are 
combined  to  form  a  new  JTRCH  (JN)  with  the  "non- 
common1'  elements  in  the  middle.   All  variations  of 
Jl  and  J2  are  deleted. 

K:  2615,   Jl:  36,   J2 :  57,  JN:  362157 
Variation  of  JN:  361257.   If  the  common  elements 
are  both  in  first  position  the  order  of  Jl  or  J2 
has  to  be  reversed  before  the  combination. 

6.  contains  the  first  and  the  last  element  of  a  JTRCH 
(J) ,  then  J  is  augmented  by  the  "non-common"  ele- 
ments forming  a  new  JTRCH  (JN) .   In  the  variations 
the  first  and  the  last  element  of  J  as  well  as 
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each  of  the  "non-common"  elements  have  to  be  in  an 
end  position  (6.  is  a  special  case  of  1.). 
K:  36752,   J:  572 
JN's:  57236,   27563 
7.   contains  the  first  and  the  last  element  of  a  JTRCH 
(Jl)  and  at  least  the  first  or  last  element  of  an- 
other JTRCH  (J2) ,  then  Jl,  J2,  etc.,  are  combined. 
Each  of  their  end  elements  and  the  "non-common" 
elements  have  to  appear  in  an  end  position  of  one 
of  the  variations  of  the  new  JTRCH  (JN) .   All  va- 
riations of  Jl,  J2,  etc.,  are  deleted. 
K:  73425,   Jl :  32,   J2 :  78 
JN's:  324578,   324587,   432785. 
The  largest  KCHAIN  contains  all  elements  of  the  trunk,  be- 
cause they  correspond  to  the  trunk  branches  present  in  the 
main  loop  (cf.  sect.  4.2),  we  will  always  be  able  to  gene- 
rate the  complete  trunk (s)  by  this  combination  procedure. 

5.4   Connecting  Main  Branches  and  Limbs 

When  the  trunk (s)  has  been  established,  the  main 
branches  and  limbs  are  connected.   In  section  4.4  this 
step  of  the  realization  procedure  is  performed  in  tabular 
form.   A  similar  "table"  is  used  in  the  computer  program, 
each  main  branch 

1.  the  number  of  loops  in  which  it  appears  is  stored 
in  the  array  MSER, 

2.  For  each  of  those  loops  the  trunk  branches  are 
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stored  in  the  array  MBRT  and  the  other  main 

branches  in  MBRB. 
Example:   Tree  branch  3  (Fig.  34)  is  the  third  main 
branch 

MSER  3=2 

MBRT  (3,1,1)  =5    MBRB  (3,1,1)  =  2   (loop  1) 

MBRT  (3,2,1)  =7    MBRB  (3,2,1)  =  2   (loop  2) 
Unique  connections  are  found  by  inspecting  the  above  table 
whether  the  conditions  of  Theorems  7a  and  7b  and  Corollary 
7  are  met  by  any  of  the  main  branches.   For  this  example 
main  branches  3  and  1  are  both  uniquely  connected  to  the 
node  between  trunk  branches  5  and  7,  either  main  branch  3 
or  1  may  be  connected  first.   Main  branch  2  is  uniquely 
connected  to  the  outer  node  of  trunk  branch  6. 
Once  all  unique  connections  have  been  made,  a  search  is  per- 
formed with  the  remaining  main  branches,  whether  they  can 
be  attached  to  the  already  connected  main  branches,  accor- 
ding to  Theorem  7b.   When  this  step  has  been  performed  and 
there  are  still  unconnected  main  branches  left,  they  are 
attached  to  trunk  or  already  connected  main  branches  such 
as  to  satisfy  Properties  1  or  3 .   Limbs  are  simply  connect- 
ed to  the  outermost  main  branch,  or  leftmost  trunk  branch 
if  there  are  no  main  branches  in  that  specific  loop. 

5 . 5   Connecting  Links 

The  links  connect  the  end  nodes  of  the  linear  paths 
formed  by  the  tree  branches  in  each  loop.   In  sect.  4.4.6 
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these  end  nodes  are  found  by  an  inspection  of  the  already- 
synthesized  tree,  represented  by  a  graph,  and  the  rows  of 
the  loop  matrix.   In  the  computer  program  an  efficient 
search  for  the  end  nodes  of  the  linear  paths  is  achieved  by 
identifying  tree  branch  patterns  for  each  loop.   These 
patterns  are  demonstrated  in  Fig.  35,  where  the  capital  let- 
ters indicate  whether  trunk  branches  (T) ,  and/or  main  bran- 
ches (M) ,  and/or  limbs  (L)  are  present  in  a  loop.   The  po- 
sition of  the  link  (dotted  line)  is  indicated  as  specified 
by  the  particular  pattern. 

This  completes  the  programmed  version  of  the  realization 
procedure.   Next,  two  typical  computer  solutions  are  given. 
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Figure  34 
Graph  corresponding  to  loop  matrix  in  Sect.  5.2 
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Figure  35 

Patterns  of  link  positions  with  respect  to  trunk  branches 
(T) ,  main  branches  (m) ,    and  links  (L) . 
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5.6   Two  Computer  Solutions 

Computer  solution  of  the  example  in  Section  4.4 
The  graph  is  shown  in  Figure  29. 
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A  general  problem.   The  graph  is  shown  in  Figure  36. 
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EVALUATION  OF  THE  TOTAL  NUMBER  OF  DIFFERENT  GRAPHS 
HAVING  AN  IDENTICAL  LOOP  MATRIX 

6. 1  Introduction 

Up  to  now  this  realization  procedure  has  been  used  to 
find  a  graph  associated  with  a  given  loop  set  matrix.   How- 
ever, there  is  another  application.   It  has  been  pointed 
out  in  section  3.1  that  there  may  exist  no  graph,  exactly 
one  graph,  or  many  graphs  corresponding  to  a  given  loop  set 
matrix.   This  realization  technique  gives  a  way  of  calcula- 
ting the  total  number  of  different  graphs  having  the  same 
loop  set  matrix.   Although  this  problem  is  at  the  present 
only  of  theoretical  importance,  it  demonstrates  the  general- 
ity of  the  realization  procedure  presented  in  this  thesis. 

6. 2  The  Concept  of  Different  Graphs 

Given  two  graphs  G,  and  G~  having  the  same  number  of 
nodes  and  labeled  edges.   Moreover,  they  shall  have  the 
same  number  of  loops.   Each  loop  in  G,  can  be  identified 
with  a  loop  in  G~,  i.e.  corresponding  loops  are  made  up  of 
the  same  edges „ 

Definition:   Two  graphs  G,  and  G~  having  the  same  number  of 
nodes,  labeled  edges,  and  identical  loops  are 
called  different  if  there  is  at  least  one  node 
which  is  not  incident  with  the  same  edges  in 
both  graphs. 
The  two  graphs  in  Fig.  37  are  different.   The  concept  of 
different  graphs  could  be  useful  in  integrated  circuit  lay- 
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outs.   If  in  Fig.  37  edges  5  and  2  are  capacitive  and 
edges  1  and  4  resistive,  the  circuit  arrangement  in  Fig. 
37b  could  be  implemented  easier  than  that  in  Fig.  37a.   To 
summarize:   different  graphs  are  physically  different,  i.e, 
the  arrangement  of  their  labeled  edges  relative  to  each 
other  is  not  the  same.   However,  these  graphs  all  have  the 
same  topological  structure. 

6. 3   Evaluation  of  the  Total  Number  of  Different  Trunks 

Before  calculating  the  number  of  different  graphs, the 
number  of  alternate  solution  for  each  realization  step  is 
determined.   Assuming  the  trunk  is  made  up  of  T  disjoint 
chains  of  trunk  branches,  i.e.  which  can  be  permuted  with- 
out violating  any  common  node  requirement  (CNR) .   Let  the 


ith  chain  be  composed  of  t.  edges  which  can  be  permuted 
within  the  chain  not  violating  any  CNR.   Then  the  follow- 
ing theorem  is  true. 
Theorem  8 :   If  the  trunk  is  composed  of  T  disjoint  chains 

and  each  chain  has  t.  permu table  edges,  then 

there  exist 

Nfc  =  (T!/2)|   I  fc  ' 

i  =1 

different  trunks. 
Proof:   This  expression  is  verified  using  the  well  known 
formulae  of  permutations.   The  factor  one  half  is  introduced 
since  two  trunks  are  not  considered  to  be  different  if  they 
have  the  same  sequence  of  branches  either  from  "left  to 
right"  or  "right  to  left". 
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Referring  to  Fig.  38  the  three  chains  can  be  permuted  in 
31/2  ways.   The  branches  within  each  chain  can  be  permuted 
in  3!  ways,  yielding  a  total  of  648  different  trunks. 

When  there  are  CNR  within  the  chains,  i.e.  some  bran- 
ches form  subchains,  the  expression  in  Theorem  8  has  to  be 
modified. 

Theorem  9 :   If  the  trunk  consists  of  T  disjoint  chains  and 
the  branches  in  the  ith  chain  form  G,  success- 
ively smaller  disjoint  groups  of  subchains, 
each  of  which  having  g. .  disjoint  components, 
then  there  exist 

T   _L 
Nt  =  (TI/2)  FT    M    g±j: 

i=l j=l 

Proof:   This  theorem  is  verified  applying  combinatorial 

mathematics.   In  Fig.  39  chains  1  and  2  can  be  permuted  in 

(2!/2)  =  1  way.   Chain  1  does  not  have  any  subchains,  thus 

there  is  only  one  CNR,  i.e.   G-.  =  1.   The  edges  1  and  2  can 

be  permuted  in  2'  ways.   Chain   2  consists  os  two  subchains, 

thus  there  are  3  CNR  within  chain  2,  i.e.   G~  =  3.   The 
3  2 

product  I    g~  .  J  is  composed  of  the  terms 

j=l    : 

g2-i  =  11      for  the  elements  in  subchain  2a 
g»~  =  3  J   for  the  elements  in  subchain  2b 
g^_  =  21       for  the  subchains  2a  and  2b  in  chain  2 
Thus  there  exist  a  total  of  (21/2)     [21 ]  [11312!]  =  24  dif- 
ferent trunks.   This  completes  the  calculation  of  the  num- 
ber of  different  trunks. 
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6.4   Evaluation  of  the  Total  Number  of  Different  Main 

Branch  Connections 

In  the  process  of  connecting  main  branches  to  the 
trunk  or  already  connected  main  branches,  there  can  arise 
four  possibilities: 

1.  One  main  branch  is  uniquely  connected  to  a  node. 

2.  Two  or  more  main  branches  are  uniquely  connected 
to  the  same  node. 

3.  One  main  branch  can  be  connected  to  either  end  of 

a  linear  path  of  trunk  branches  and/or  already  con- 
nected main  branches. 

4.  More  than  one  main  branch  can  be  connected  as  in  3, 
Clearly  case  1  does  not  yield  alternate  solutions. 

Assuming  there  are  M  groups  of  main  branches  uniquely  con- 
nected to  the  same  node,  the  ith  group  having  m.  elements. 
Assuming,  moreover,  there  are  N  groups  of  main  branches 
that  could  be  connected  to  either  one  of  the  two  end  nodes 
of  a  linear  path,  the  jth  group  having  n.  elements.   Then 
following  theorem  holds: 

Theorem  10:   If  the  main  branches  can  be  uniquely  connected 

in  M  groups  to  the  same  node,  such  that  each 

elements  of  the  group  may  be  connected  first, 

and  N  groups  of  main  branches  can  be  connected 

to  either  of  two  end  nodes  of  a  linear  path, 

then  there  are 

M        N 


N   =     m. !      (n.+l) : 
m    i-1   x    j=l   3 
different  ways  to  connect  these  main  branches. 
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Proof:   Referring  to  Fig.  40a  the  3  main  branches  can  be 
connected  in  31    different  ways  to  node  v.   From  Fig.  43b  it 
follows  that  main  branch  4  can  be  connected  in  two  ways  to 
end  nodes  a  or  b.   Fig.  40c  reveals  that  there  are  4x3' 
different  ways  to  connect  the  main  branches  5,6,  and  7  to 
either  one  of  the  end  nodes  c  and  d.   If  all  main  branches 
are  connected  to  one  node,  3!  different  connections  result. 
Alternately,  one  branch  can  be  selected  in  3  ways  and  then 
connected  to  one  node.   The  remaining  2  main  branches  are 
then  attached  in  21    different  ways  to  the  other  node. 
This  again,  yields  31    different  possibilities.   This  com- 
pletes the  calculation  of  the  different  possibilities  to 
connect  the  main  branches. 

6. 5   Evaluation  of  the  Total  Number  of  Different  Limb  and 

Link  Connections 

Finally  the  number  of  alternate  connections  for  limbs 
and  links  is  given  by 

Theorem  11:   If  there  are  L  loops,  the  ith  loop  having  I. 

limbs,  then  the  limbs  and  the  links  can  be  con- 
nected in 

L 


Ni  =  I       (t±+i): 

i=l 

Proof:   Referring  to  Fig.  41  there  are  31   ways  to  connect 
the  three  limbs,  assuming  the  link  is  connected  last.   How- 
ever, the  link  can  be  placed  in  four  different  positions, 
thus  there  are  a  total  of  4x3  I   different  ways  to  con- 
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nect  the  three  limbs  and  the  link.   This  completes  the  num- 
ber of  different  ways  limbs  and  links  can  be  connected. 

6.6   Evaluation  of  the  Total  Number  of  Different  Graphs 

By  the  definition  of  the  term  "different  graph"  each 
of  the  alternate  connections  yields  a  different  graph.  Thus 
the  following  theorem  holds: 

Theorem  12 ;   If  the  realization  procedure  presented  reveals 
that 

1.  The  trunk  consists  of  T  disjoint  chains, 
the  ith  chain  forming  G.  successively 
smaller  disjoint  groups  of  subchains  each 
having  g. .  disjoint  components, 

2.  the  main  branches  can  be  uniquely  connected 
in  M  groups  to  one  node  and  attached  in  N 
groups  to  either  one  of  two  end  nodes, 

3.  there  are  L  loops,  the  rth  loop  having  I 
limbs  then  there  exist 


N  = 


T    i  N 


<T:/2>  l.j  H  'iJ^U  'V  I   l<V1):^l  K+1>;1 
1=1  3=1    J    k=l    p  =1  v  r=l 


different  graphs  corresponding  to  a  given 
fundamental  loop  or  cut  set  matrix. 
The  unknowns  present  in  this  expression  can  be  obtained 

from  the  computer  program.   The  common  node  requirements  are 

listed  in  the  output  for  this  purpose. 
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CONCLUSION  AND  SUGGESTION  FOR  FURTHER  RESEARCH 
The  efficiency  and  generality  of  a  new  realization 
procedure  for  a  fundamental  loop  or  cut  set  matrix  has 
been  demonstrated.   The  concept  of  trunk  branches,  main 
branches,  and  limbs  proves  to  be  useful  in  describing 
graphs.   The  idea  of  physically  different  networks  is 
introduced  which  should  prove  to  be  useful  in  finding  or 
classifying  equivalent  networks. 
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APPENDIX  I 

Glossary  of  Topological  Terms 

The  definition  of  topological  terms  used  in  this  thesis 

are  summarized  below. 

Adjacent:  Two  edges  are  adjacent,  if  they  have  a  vertex  in 
common. 

Basic:   See  Fundamental 

Branch:  See  Tree  Branch 

Connected:  A  graph  is  connected  if  there  exists  a  path  be- 
tween any  two  nodes;  a  connected  graph  cannot  contain 
any  isolated  node. 

Co-tree:   The  edges  of  the  graph  which  are  not  part  of  the 
tree  form  the  complement  of  the  tree  or  co-tree. 

Cut  Set:   A  set  of  edges  forms  a  cut  set  if  the  removal  of 
this  set  "cuts"  the  graph  into  separate  parts,  however, 
removing  all  but  one  of  these  edges  leaves  the  graph 
connected. 

Edge:   Two  distinct  points  (end  points)  and  the  line  seg- 
ment joining  them  form  an  edge. 

Fundamental  Loop:   A  fundamental  loop  consists  of  one  link 
and  the  tree  branches  connecting  its  nodes. 

Fundamental  Cut  Set:   A  cut  set  which  besides  links  cuts 

only  one  tree  branch. 

Fundamental  Loop  Set  Matrix 

links  tree  branches 


B 


U       F  loops 

b. .  =  1  if  edge  j  is  present  in  loop  i. 
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Fundamental  Cut  Set  Matrix 

links  tree  branches 
Q    =     [_-F        U  cut  sets 

q.  .  =  1  if  edge  j  is  cut  by  the  cut  i. 

Graph:   A  "linear  graph"  or  just  graph  is  a  collection  of 

edges  with  no  single-edge  loops. 
Incident:   A  node  and  an  edge  are  incident  if  the  node  is 

an  endpoint  of  the  edge. 

Incidence  Matrix 2 

edges 

-| 

nodes 


A   =     a. 

a. .=  1  if  edge  j  and  node  i  are  incident 
ID 

Incidence  Set :   The  set  of  edges  incident  at  one  node  from 
an  incidence  set. 

Isomorphisms  :   If  a  graph  G-.  can  be  broken  at  a  single  node 
into  two  disjoint  subgraphs  or  two  disjoint  subgraphs 
are  joined  at  a  single  node,  then  the  resulting  graph 
and  G,  are  1-isomorphic .   If  a  graph  G,  is  cut  at  two 
nodes  into  two  disjoint  subgraphs,  one  of  them  is 
turned  around  and  then  joined  with  the  other  at  the 
two  nodes,  G,  and  the  resulting  graphs  are  2-isomor- 
phic. 

Limbs   A  tree  branch  which  appears  in  only  one  loop. 

Linear  Graph;   See  Graph. 

Link;   An  edge  of  a  graph  which  is  not  a  tree  branch,  is  a 
link.   The  links  form  the  co-tree. 
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Loop:   A  connected  subgraph  in  which  exactly  two  edges  are 

incident  with  each  node. 
Main  Branch:   A  tree  branch  which  appears  in  more  than  one 

loop  and  is  not  part  of  the  trunk. 
Main  Loop:   The  loop  of  a  graph  containing  the  maximum 

number  of  those  tree  branches  which  are  common  to  more 

than  one  loop. 
Node:   The  endpoint  of  an  edge  is  called  a  node  or  vertex. 

Isolated  points  of  a  graph  are  called  "isolated  nodes". 
Oriented:   If  an  edge  is  assigned  a  direction,  the  edge  is 

oriented.   If  a  graph  is  made  up  of  oriented  edges, 

the  graph  is  called  oriented. 
Path:   A  linear  path  or  just  path,  between  nodes  i  and  j  is 

a  sequence  of  edges.   Each  node  of  the  path  is  shared 

by  exactly  two  edges  except  nodes  i  and  j  which  form 

the  endpoints  of  the  path. 
Separable  Graph:   A  connected  graph  is  separable  if  it  con- 
tains at  least  one  subgraph  which  has  only  one  node  in 

common  with  its  complement. 
Single-Edge  Loop:   Coalescing  the  two  nodes  of  an  edge 

generates  a  single-edge  loop. 
Tree:   A  connected  subgraph  containing  all  nodes  of  a  graph 

but  no  loops  is  called  a  tree  of  a  graph. 
Tree  branches:   The  edges  of  a  graph  which  are  part  of  the 

tree  are  called  tree  branches  or  just  branches. 
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Trunk:   The  linear  path  which  is  formed  by  the  tree  branches 
appearing  in  the  main  loop  (omitting  the  limbs)  is 
called  the  trunk  of  the  tree. 

Trunk  Branches :   The  tree  branches  which  form  the  trunk  are 
called  trunk  branches. 

Unique  Connection:   If  the  position  of  a  main  branch  rela- 
tive to  the  trunk  and/or  connected  main  branch  is  ex- 
actly specified,  the  main  branch  is  uniquely  connected 
with  respect  to  those  branches. 

Vertex:   See  Node. 
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APPENDIX  II 

Transformation  of  Current  and  Voltage  Sources 

For  each  electric  network  there  exists  an  associated 
linear  graph.   This  graph  can  be  found  by 

1.  open  circuiting  ideal  current  sources, 

2.  short  circuiting  ideal  voltage  sources, 

3.  and  replacing  passive  elements  by  edges. 

If  there  are  branches  in  the  network  which  consist  of  ideal 
sources  only,  these  sources  can  be  transformed  such  that 
each  of  the  ideal  current  sources  is  parallel  to  a  passive 
element,  and  each  of  the  ideal  voltage  sources  is  in  se- 
ries with  a  passive  element.   This  transformation  techni- 
que is  taken  from  Reference  4.   It  is  assumed  that  there 
are  no  loops  formed  by  independent  voltage  sources.   Assum- 
ing this  were  not  the  case,  one  of  these  sources  had  to  be 
dependent  by  KVL,  which  contradicts  the  assumption. 

Since  voltage  sources  do  not  form  loops,  an  edge  con- 
sisting of  a  voltage  alone  can  be  identified  as  a  tree 
branch  (Fig.  42a) .   Removal  of  this  tree  branch  cuts  the 
tree  into  two  separate  parts.   The  voltage  source  is  placed 
in  series  with  each  element  corresponding  to  the  tree  bran- 
ches in  one  of  the  separate  parts.   The  two  nodes  of  the 
removed  edge  are  then  coalesced  (Fig.  42b) .   This  transfor- 
mation does  not  change  the  network  equations.   This  is 
verified  by  writing  KVL  for  the  loops  containing  the  tree 
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branches  into  which  the  voltage  sources  have  been  inserted, 
or  KCL  for  the  two  nodes  which  have  been  coalesced. 

Since  current  sources  do  not  form  cut  sets,  and  edge 
consisting  of  a  current  source  alone  can  be  identified  as  a 
link  (Fig.  43a) .   The  current  source  is  placed  in  parallel 
with  each  element  corresponding  to  those  tree  branches  which 
form  a  loop  with  that  link  (Fig.  43b) .   This  transformation 
does  not  change  the  network  equations.   This  is  verified  by 
writing  KCL  for  the  two  nodes  of  that  link,  or  KVL  for  the 
voltage  between  the  two  nodes  of  that  link. 

Thus  each  branch  of  an  electric  network  can  be  repre- 
sented in  its  most  general  form  as  in  Fig.  44,  where  the 
sources  (v  ,  and  j  ,  )  may  be  absent,  the  passive  element  R, 
is  always  present. 
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Figure  42 
Voltage  source  transformation 


Figure  43 
Current  source  transformation 
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Figure  44 
Typical  branch  of  an  electric  network 
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APPENDIX  III 
Computer  Program 

C  THIS  PROGRAM  REALIZES  A  FUNDAMENTAL  LOOP  SET  MATRIX, 
C   INPUT: 
C      1  )  ROWS  AND 
C      2)  COLUMNS  OF 

C         THE  AUGMENTFD  F-MATR I X , WHOSE  ELEMENTS  IN  COLUMN  1 
C         CORRESPOND  TO  LINKS, AND  WHOSE  ELEMENTS  IN  ROW  I 
C         CORRESPOND  TO  TREE  BRANCHES. 
C      3)  THE  AUGMENTED  F-MATPIX,ONE  ROW  PER  CARD. 
C  OUTPUT: 

C       1)  AUMGNENTED  F-MATRIX  PUT  IN, 
C       2)  PARTITIONED  F-MATRIX, 
C       3)  COMMON  NODE  REQUIREMENTS, 

C       4)  NUMBER  OF  DIFFERENT  TRUNKS  WHICH  ARE  ARE  IMPORTANT 
C  FOR  THE  ALGORITHM  USED  IN  THIS  PROGRAM, 

C       5)  ONE  GRAPH. IF  ALTERNATE, HOWE VER  EQUIVALENT, 
C  REALIZATIONS  EXIST, SET  THE  VARIABLE  "NOGRAF"  TO 

C  TO  THE  DESIRED  NUMBER. IF  THE  PROGRAM  RECOGNIZES 

C  FEWER  ALTERNATIVES. ONLY  THESE  ARE  PRINTED. 

IMPLICIT  INTEGER  (A-Z) 

DIMENSION    F (30.30) , LIMB(2C) ,NCTLIM(20) ,MAINBR(20) 

DIMENSION    TRUNBR(20) 

DIMENSION    MEND(30) , GRAPH ( 60 , 3 ), TRUNK ( 20 , 20 ), FpART ( 3C ) 

DIMENSION  MSER(30),MBR9(3C,10,1C  )  ,MBRT(3C  ,10,]"!) 

DIMENSION  CONECK40) 

DIMFNSION  MTS(5,4) , MT ( 5 ) , MBGS ( 5 ) , MAMBIG( 9 , 9 ) 

DIMENSION  MSLET(30,1C  > , MSLE8( 30 , 1C ) ,MCOM(5,5) ,MC<9) 

DIMENSION  JTR CH( 30, 3C), LENGTH (2T ) , JLAST ( 10  I , J F IRST <  IC  > 

DIMENSION  JTEMP1 (30) 

DIMENSION  JADD(30,9).MUNIQ<9,4) ,MCOMP(10)  ,JINTLE(I) 

DIMENSION  KSTORE  ( 2C I .KCHAI N( 20 , 10 , 20 ) , KC ( 10 ) 

DIMENSION  JTEMP(30) .JSEQ(30) 

DIMENSION  KTRCH(10,1C ),KLE(2C ), JSTOR( 5,10) 

COMMON  /CH/  KCHAIN, KSTORE, KC 

COMMON/PT/JTRCH, LENGTH, JTWO, JLAST, JF IRST, J ADD, 
1JTEMP1, JT,JF, JL 

COMMON/ LLO/  ROW,COL,F,TR,TRUNBR,Me,TRUNK,MBR,MAINBR, 
1MEND,MGR,GPAPH,L,LIMB,ME 

COMMON  /NOD/JDELET(40), JD,COMB(^0),C,J 

READ  (5,299)  ROW 

READ  (5.299)  COL 
2^9   FORMAT  (12) 

DO  297  IN*I, ROW 

READ  (5,29P>  ( F ( I N, I ) , I =1 ,COL ) 

297  CONTINUE 

298  FORMAT  (3012) 
NOGRAF=l 

C  PARTITION  F  MATRIX 
C  FIND  LIMBS 

L=0 

NL=0 

00    2       11=2, COL 

COLSUM=0 

DO    3    12=2, ROW 

3  COLSUM  =  COLSUM+F(  12,  ID 
IF  (COLSUM.GT.l)GOTO  4 
L=L  +  1 

LIMB(L)=I1 
GOTO    2 

4  NL  =  NL-»-l 
NOTLIM(NL)  =  U 

2  CONTINUE 

C    FIND    MAINLOOP 
ITEST=0 

DO    5       I3=2,ROW 
ROWSUM=0 
DO    1       14=1, NL 
C0LF=N0TLIM(I4) 
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5 
C  FT 
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3.1 


3.)  5 


3  )7 
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C  FI 


10 


9 
C  Fl 

C   B 


ROW 
IF( 

TRO 
ITE 
CON 
TR  = 
MBR 
NO  TR 
DO 
COL 
IF 
MBR 
MAI 
GOT 

TR  = 

TPU 

CON 

WRI 

DP 

WRI 

CON 

FPA 

IF( 

no 

FPA 
CON 
IF( 
DO 
FPA 
CON 
IF( 
DO 
FPA 
CON 
WRI 
DO 
WPI 
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KSU 
DO 
COL 
IF 
KSU 
KST 
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Y  COM 
IF 
DO 
M*  I 

MFL 
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KT  = 
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ROWSUM.IE.I TEST) 

W=I? 

ST=ROWCHM 

TINUE 


COLF) 
GOTO  5 


=  0 

UNK  BP 
6   15  = 
F=NOTL 
(F(TRO 
=  MBR*-1 
NRR< MB 
0  6 
TR*! 
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TINUE 
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W,COLF).EQ.l )  GOTO  7 

R)=COLF 
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5f  ) 

=  1 ,ROW 
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r)GO  TO  3  0 
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5f )=TRUNBR 
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1=1  ,  MBP 

R  +  N51 ) =MAI 
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)GOTO  1C 
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DO  13   K4=2,R0W 

IF(F(K4,MAINF).EQ.C)G0T0  13 

KSUM=0 

DO  14   K5=1.TR 

COLF=TRUNBP(K5) 

IF(F(K4,C0LF>.E0.0)GOTO  14 

KSUM=KSUM+1 

KSTORE(KSUM)=FU,COLF) 

14  CONTINUE 

IF  (KSUM  -1)  13,15,17 

15  MFLAG=MFLAG  «-l 

IF  (MFLAG.GT.J.  )GOTO  16 
MTEMP=KSTOPE(l) 
GOTO  13 

16  MTEMP1=KSTCRE(JJ 

IF  (MTEMP.EO.MTEMP1)  GOTO  13 
KSUM=2 

KSTORE(l)=MTEMP 
KSTORE(2)=MTEMPl 
CALL  CHAINS  (KSUM) 
MFLAG=0 
MFLAG2=1 

17  KT=KT*1 

DO  18  K6=1,KSUM 

18  KTRCH(KT,K6)=KST0RE(K6» 
KLE(KT)=KSUM 

13    CONTINUF 

IF  (MFLAG.GE.1.AND.MFLAG2.EQ.C )     GOTO  19 
GOTO  2^ 

19  KT  =  KT«-1 
KTRCH(KT,1)=MTEMP 
KLE(KT)=1 

20  CONTINUE 

IF  (KT.LE.l)  GOTO  12 

KCHFCK=KT-1 

DO  21  K7=1,KCHECK 

K8=KLE(K7) 

K9=K7+l 

DO  22  K10=K9,KT 

KSUM=0 

K11=KLE(K1C» 

DO  23  K12=1,K1! 

KNOT=0 

DO  24  K13=1,K8 

IF  (KTRCH(K10,K12).E0.KTRCH(K7,K13) )  GOTO  24 

KN0T=KN0T+1 

24  CONTINUE 

IF    (KN0T.NF.K8)    GOTO    23 
KSUM=KSUM+1 

KSTORE  (  KSUM  )=KTRCH(  KICK  12) 
23         CONTINUE 

IF    (KSUM.EQ.C)     GOTO    22 

CONTINUE 

DO    25    K14=1,K8 

KSUM=KSUM+1 

25  KST0RE(KSUM)=KTRCH(K7,K14) 
CALL    CHAINS(KSUM) 

22    CONTINUE 

21  CONTINUE 
12  CONTINUE 
11    CONTINUE 

WRITE(6,357» 

DO  310  N59=2,TR 

IF(KC(N59).EQ.C)G0T0  310 

N60=KC(N59) 

DO  311  N61=1,N60 

WRITE (6,3091  ( KCHAIN( N59 ,N6 1 , I ) ,I=1,N59) 
311   CONTINUE 
310   CONTINUE 

JTWO=0 

J=0 
COMBINE    TRUNKCHAINS    TO    LARGER    CHAINS 
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CALLED  JTRCH 

DO  30  J1=2,TR 

IF  (KCUl).EQ.C)  GOTO  30 

JKC=KC( Jl) 

DO  31  J2=1,JKC 

IF  (J.NE.O)  GOTO  32 
^3    JKL  =  1 

DO  34  J3=1?J1 
34    JTEMP1  <J3)=KCHAIN( J]  ,J2fJ?> 

JT  =  J1 

GOTO  35 

32  CONTINUE 

IF  U.EO.C)  GOTO  610 

DO  608  J6C=1.J 

J6l=IABS(LENGTH( J60)) 

LENGTH( J6C)=J61 
6  38   CONTINUE 
610   CONTINUE 

DO  29  JCN  =  ltJl 
29    JTEMPK  JCN)=KCHAIN(J1  ,J2,JCN) 

JT  =  J1 

JD  =  ^ 

CALL    NODUP 

JT  =  0 

JINT=0 

JTWO=0 

JFLAG1=0 

JFLAG2=0 

JMOP=0 

JNOTFL=0 

JF=0 

JCONT=0 

JL=0 
CHFCK    WHICH    JTPCH'S    CAN    BE    GROUPED    TOGETHER 

DO    36    JC4=1 ,C 

J4=C0MB< JC4) 

JLE  =  IABS(LFNGTHU4) ) 

LENGTH! J4)=JLE 

JF1=0 

JLI=0 

JCOUNT=0 

DO  37  J5=1,J1 
37    JTEMPC J5)=KCHAIN( Ji,J2,J5) 

JN=0 

JFLAG=r 

JSTART=1 

DO  38  JV=1,J1 

IF( JTEMP( jv).EO. JTRCH (J4,l) I  JF1=JV 

CONTINUE 

IF< JTEMP(JV). EQ. JTRCH <J4,JLE> )  JL1=JV 
3  3    CONTINUE 

IF  <JF1#E0.0.  0R.JL1.EQ.C)  GOTO  26 

LENGTH! J4)=-JLE 

JM0R=JM0R+1 
26    CONTINUE 

IF  (JF1.NE.0)  GOTO  39 

CONTINUE 

IF  (JL1.NE.0)  GOTO  40 

GOTO  41 

40  JL=JL+1 
JFLAG=1 
JLAST1 JL)=J4 
JTEMP( JL1)=JTEMP(1) 
GOTO  4? 

33  JF  =  JF-H 
JFLAG=2 
JFIRST( JF)=J4 
JTEMPC JF1)=JTEMP(1) 

42  CONTINUE 
JC0UNT=1 
JSTAPT=2 

41  CONTINUE 
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JST=JSTAPT 

00  43  J6=JST,JLE 

JR  =  J6 

IF( JFLAG.EO. 1)  JR=JLE+1-J6 

CONTINUE 

00  44  J7=JSTART,J1 

IF( JTE MP (J7).EQ. JTRCH (J4,JR» )GOTO  4  5 

44  CONTINUE 
JN=JN*1 

IF( JN.FO. ( JLE-JST) )  JNOTF L= J NOT FL  +  1 
GOTO  43 

45  JTEMP( J7)=JTEMP(JSTART> 
JSTART=JST£RT-»-? 

IF( JCOUNT.FQ.OJGOTO  46 

JC0UNT=JC0UNT+1 

IF( JC0UNT.FQ.J6)G0T0    43 

JD=JD*i 

JDELET( JD)=J4 

GOTO  36 

46  JC0UNT=J6 
43    CONTINUE 

IF(  JSTART.FO.  (  JU1II  GOTO  47 

CONTINUE 

IF  (JFLAG.Nt.O)  GOTO  48 

IF  ( JNOTFL.GE. J>  JFL&G2=1 

GOTO  36 

48  JINT=JINT+! 
JIC=0 

DO  49  J8  =JSTART,J1 
JIC=JIC+1 

49  JSTOR( JINT, JIC)=JTEMP(J8> 
JINTLE< JINT)=JIC 

JSEO( JINT)=J4 
GDTO  36 

47  JCONT=l 

IF  (JFLAG.EO. 21  GO  TO  28 

CONTINUE 

IF( JFLAG.EO. II  GO  TO  2? 

CONTINUE 

IF( JFLAG.EO. 01  GO  TO  36 
27    JL=JL-1 

GOTO  36 
23  JF=JF-l 
36    CONTINUE 

IF( JD.EQ.OIGOTC  51 

50  CONTINUE 
JNEW=0 

DO  52  J9=1,J 

TEST=J9 

NEW=JNEW 

CALL  JTRONE  ( NEW, TEST , 652 > 

IF( JD.EO.CIGOTO  52 

DO  53  J10=1,JD 

IF  (  J9.E0.  JDELET(  J10HGOTO  52 

53  CONTINUE 
JNEW=JNEW+1 
J!1=IABS(LENGTH( J9) ) 
DO  54  J12=1,J11 

54  JTRCH(JNEW,J12)=JTRCH( J9t J12) 
LENGTH( JNEW)=J11 

52         CONTINUE 
J=JNEW 
IF< JFLAG1.E0.1IG0T0  31 

51  CONTINUE 

C  KCHAIN  IS  COMPLETELY  CONTAINED  IN  A  JTRCH 

IF  (JCONT.EQ.l)  GOTO  31 

CONTINUE 

IF  ( JFLAG2.EQ.1)  GOTO  33 

CONTINUE 
C  THREE  OR  MORE  JTRCH  ARE  COMBINED 

IF( (JF+JL+JM0R).LE.2)G0  TO  55 

KTOTLE=0 
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DO  56   K3C  =  1 .JINT 
K31=JSEQ(K30) 
K32=IABS(LENGTH(K31> ) 

56  KTOTLE=KTOTLE+K22 
JNTR=0 

IF  (KTOTLE.EO. Jl)  GOTO  57 

J  =  J*1 

DO  58   K33  =  1,J1. 

KNOT=0 

DC  59   K34=1,JINT 

K35  =  JSFO(K?M 

K36=IABS(LFNGTH(K35) ) 

DO  60   K37=1,K36 

IF  (KCHAIM Jl , J2,K33).EQ. JTRCH< K35 , K37 ))  GOTO  60 

KN0T=KN0T+1 
-O    CONTINUE 
59    CONTINUE 

IF  (KNOT.NE.KTOTLE)  GOTO  58 

JNTR=JNTR+1 

JTRCH < J,JNT9|=KCHAIN( Jl ,J2,K33> 
58    CONTINUE 

JINT=JINT+1 

JSEQ( JINT)=J 

LENGTHl J)=JNTR 

57  CONTINUE 

DO    61       K38=1,JINT 

K39=JSE0(K38I 

K40=IABS(LENGTH(K39)  ) 

JD=JD*J 

JDELET( JD)=K3^ 

J  =  J+1 

LONG=0 

DO  62   K41=l»K4r 

LONG=LONG+l 

62  JTRCH( J,L0NG)=JTRCH(K39,K4l ) 
K42=JINT-1 

DO  63   K43=lfK4? 

KL=K38*K43 

IF  (KL.GT.JINT)  KL=KL-JINT 

K44=JSEQ(KL> 

K45=IABS(LENGTH(K44) ) 

DO  63   K46=1,K*5 

LONG=LONG+l 

63  JTRCH( J,L0NG)=JTRCH(K44,K46» 
LENGTH( J)=LONG 

61    CONTINUE 

IF  (JNTR.LF.2)  GOTO  64 

JPREV=J 

JPLE=LENGTH( J) 

JMOREP=( JNTR-2) /2 

JODD=C 

I^(  (  (JNTR-T  )/2).EQ.  (JNTR/2)  )  JODD=JNTR- JNTP /2 

CONTINUE 

IF  (JODD.NF.C)  JMOREP=JMOPEP+! 

CONTINUE 

DO  65   K47=l , JMOREP 

J=J  +  1 

K49=JSE0( JINT) 

K50=LENGTH(K49) 

K51  =  JTRCH(K49,  UK47) 

K52=JTRCH(K49, K50-K47) 

IF  (K47.E0. JMOREP. AND, JODD.NE.O)  K5l= JTPCH( K4Q, JCDD I 

CONTINUE 

IF  (K47.EQ. JMOPFP.AND. JODD.NE.H )  K52= JTRCH (K49, 1 ) 

JTRCH< J,1)=K51 

COUNT=l 

DO  66   K53=l ,K5C 

IF  < JTRCH (K49,K5  3).E0. K51. OR. JTRCH ( K49 , K53 ) . EQ.K 52) 
1  GOTO  66 

C0UNT=C0UNT+1 

JTRCH( J,C0UNT)=JTRCH(K49,K5  3) 
66    CONTINUE 
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K54=K50+1 

DO  67   K55=K54,JPLE 

COUNT=COUNT+l 

67  JTRCH( J,COUNT)=JTRCH( JPREV,K55) 
JTRCH( J,JPLE>=K5^ 

6S    LENGTH( J)=J°LE 
64    JFLAGl=l 

GOTO  5: 
C  TWO  JTRCH  ARE  COMBINED, JKL  INDICATE  HOW 
55    J22=JINTIE(1) 

IF  ( JINT.GT.l )  GOTQ  6B 

JT=J22 

DO  69  J23=1,J22 
69    JTEMP1(J23)=JSTCP(1 ,J23> 

GOTO  70 

68  J24=JINTLE(2) 
JT  =0 

DO  71  J25=1,J2? 

DO  72  J26  =  ?.  ,  J24 

IF  (  JSTORd,  J25).NE.JST0R(2,J26  )  )  GOTO  64 

JT=JT+! 

JTEMPK JT)=JSTnp(i ,  J^5) 

72  CONTINUE 
71  CONTINUE 
7^    CONTINUE 

IF( JNOTFL.GE.J) JKL=1 
CONTINUE 

IF( JF.EQ.H. AND. JL.EQ.ll JKL =2 
CONTINUE 

IF(  JF.EO.l.  AND.  JL.EQ.i;)  JKL  =  3 
CONTINUE 

IF( JF.EO.l. AND. JL.EO.l > JKL=4 
CONTINUE 

IF( JF.EO.O. AND. JL.EQ. 2) JKL  =  5 
CONTINUE 

IF( JF.EQ.2.AND. JL.EQ.C) JKL=6 
CONTINUE 

IF  (JKL.GF.4)  JTWO=JTWO+i 
35    CONTINUE 

CALL  PERMUT( JKL , JNLA ST, 674) 
CONTINUE 
C  A  NEW  JTRCH  IS  GENERATED  FROM  A  CHAIN 
DO  73  J27=! , JNLAST 
J  =  J  +  1 

LENGTH! J) =JT 
DO  73  J28=1,JT 

73  JTRCH< J,J2S)=JADD( J27,J28) 
GOTO  31 

74  MULT=JNLAST 

IF  ( JNLAST. EQ.O)  MULT=1 

FWD!.  =  MULT 

J29=JLAST( 1 ) 

JD=JD*1 

JDELET( JD)=J29 

IF  (JTWO.EQ.O)  GOTO  75 

J30  =  JFIRST(1.» 

JD=JD«-1 

JDELET( JD)=J3C 

IF  (LENGTH(  J3D.LT.C  )  FWD1  =  MULT=>? 

J32=IA3S(LENGTH(J30)> 

75  TIMES=1 

IF  (LENGTH* J29).LT.C)  TIMES=2 
J31=IA6S<LENGTH< J29)) 
JAD=0 
JFC=0 

DO  76  J33=! , TIMES 
DO  77  J34=l ,FWD1 
J=J*1 

DO  78  J35=l , J31 
JLD=J35 

IF  (J33.E0.2)  JLD=J31+1-J35 
73    JTRCH< J, J35)=JTRCH< J29,  JLD) 
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IF  ( JNLAST.SO.n)  GOTO  79 

JAD=JAD«-1 

IF  ( JAO.GT.JNLAST)  JAD  = JAD-JNLAST 

CONTINUE 

DO  80  J36=1,JT 
3 J    JTRCHC J,J31+J36)=JADD( J  AD, J 36) 
70    CONTINUE 

IF  (JTWO.FO.O)  GOTO  81 

JFC=JFC*1 

IF  ( JFC.GT.FWD2)  JFC= JFC-FWD2+1 

DO  82  J37=?  ,  J3? 

JFD=J37 

IF  (JFC.LE.i--l  JFD  =  J32  +  1-J37 
8?    JTRCH( J, J31+JT+J37)=JTRCH(J3C ,JFO) 
81    CONTINUE 

LENGTFM  J)=J31+JT 

IF  (JTWO.NE.O)  LENGTH! J)=LENGTH( J)  ♦  J32 
77    CONTINUE 
76    CONTINUE 

JFLAG1=1 

GOTO  50 
31    CONTINUE 
30    CONTINUE 

T=0 

DO  83  J42=1,J 

IF  (LENGTHf J42I.NE.TR)  GOTO  83 

T=T+1 

J43=LENGTH( J42) 

DO  84  J44=1,J43 
84    TRUNK(T, J44(=JTPCH( J4?,J4^) 
33    CONTINUE 

WRITE(6,358)  T 

DC  312  N62=1,T 

WRITE (6,  309)  <TRUNK(N6  2,I ),I=J ,TP ) 
312   CONTINUE 
C  POR  EACH  MAIN  BRANCH  THE  LOOPS  IT  APPEARS  IN 
C      ARE  LISTED, ALSO  THE  OTHER  TRUNK  AND 
C    MAIN  BRANCHES  PRESENT  IN  THFSE  LOOPS 

DO   35  Ml=3 ,MBR 

M2=MAINBR(M1 ) 

M  c  sp 

DO  '86  M3=2,RGW 

IF  (F(M3,M2).E0.C)  GOTO  8  6 

MS=MS+1 

MST=C 

DO   87  M4=! ,TR 

M5=TRUNBR(M4) 

IF  (F<"»3,M5).EQ.0I  GOTO  87 

MST=MST+1 

MBRT(MJ.  ,MS,MST)=F(  1  ,M5) 
^7    CONTINUE 

MSLET(M1 f MS )=MST 

MSB=0 

DO   88  M6=] ,MBR 

M7=MAINBR(M6) 

IF  (F(M3.M7).E0.C.0R.M6.E0.M1)  GOTO  33 

MSB=MSB+1 

MflRB(Ml ,MS,MSB)=F( 1,M7) 
33    CONTINUE 

MSLEB(M1,MS)=MSB 
36    CONTINUE 

MSER(M1 )=MS 
33    CONTINUE 

DO   89  M3=J ,T 

IF  (M8.GT.NGGRAF)  GOTO  89 

CONTINUE 

DO  217  MREF=1 ,MBP 
217   MSER(MREP)=IABS(MSFR(MRFP) ) 

ME  =  ^ 

MU  =  0 

MA  =  0 
C  FIND  UNIQUE  CONNECTIONS  TO  A  NODE  BETWEEN 
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TWO  TRUNK  BRANCHES 

DO   90  M].0  =  1,MBR 

MC0L=MAINBP(M10» 

MH=  MSER(MIO) 

MFIRST=0 

TRl=TR-l 

DO   91  M14=l,TPl 

M15=TRUNK(MR,M1A) 

NEX  =  TRUNMM8,M14+1) 

MFS=0 

MSS=0 

ISUNIO=0 

LWTBR=0 

DO   92  M16=1,M11 

IF  (^SLET(M!0,MI6).E0,0I  GOTO  92 

LWTBR  =  LWTBF«-1 

MFIRST=C 

MSECON=0 

M17=MSLET(M10,M16) 

DO      93    M18=lfM17 

IF  (M15.E0.MBRT<MlC,M16fM18))  MFIRST=1 

CONTINUE 

IF  (M15.E0.MBRT(M10,M16,M18) )  MFS=MFS+1 

CONTINUE 

IF  (NEX.EQ.MBRT(M10,M16,M18))  MSECPN  =  3 

CONTINUE 

IF  (NEX.E0.MBRT(M10,M16,M18) )     MSS  =  M$S«-1 
93    CONTINUE 

IF  <<MFIRST*MSECON).NE.l)   GOTO  91 

ISUNIQ=ISUNICH1 
92         CONTINUE 

IF  (MFS.EQ.O.OR.MSS.EO.C )  GOTO  91 

CONTINUE 

IF  (  ISUNIQ.EO.LWTBR.AND.LWTBR.NE.O)  GOTO  94 
91    CONTINUE 

GOTO  90 
94-    CONTINUE 

MT(1)=0 

MT(2)=0 

DO  106  M19=1.M11 

IF  (MSLET(M10,M19).EO.O)  GOTO  106 

M20=  MSLET(M10,M19) 

DO   95  M21=1,M20 

IF  <MBRT(M!0,M19f M21).EQ.M15I  GOTO  96 

CONTINUE 

IF  (MBRT(M10.M19,M21)«NE,NEX)  GOTO  95 

MT(2)=MT(2)+I 

IF  (MSLEB(M10,M19).E0.0»  GOTO  1C5 

MY=MT(2) 

MTS(2tMY)=M19 

GOTO  106 
96    MT(1)=MT(1)-H 

IF  (MSLEB(Ml0tM]9).EO,0)  GOTO  105 

MX=MT(1) 

MTS(lfMX)=M19 
95    CONTINUF 
106   CONTINUE 

MBG=0 

DO   97  M22=l,2 

NBG=0 

M23=MTS(M22.1) 

M24=MSLEB(M10,M23) 

DO   98  M25=1,M2A 

MCOT=0 

IF(MT(M22).EQ.l)  GOTO  100 

M26=MT(M22) 

DO   99  M27=1,M26 

M28=MTS(M22,M27) 

M2<'  =  MSLEB(M10,M28) 

DO   99  M30=1,M29 

IF(MBRB(M1C ,M2  3,M25).E0.MBRB(M1C, M2  8, M29) ) MCnT=MCOT+l 
99    CONTINUE 
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I  )0 


58 
57 


102 
103 
101 


104 


10  5 
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IF(MCOT.EQ. ( 
GOTO  98 
NBG  =  NB04-J 
MC0M(M??,NP0 
MC(M22)=NBG 
CONTINUE 
CONTINUE 
IF(MC(1).EQ. 
M31=MC(1) 
M32=MC(2) 
DO  1C1  M33=l 
DO  102  M34=l 
IF(MC0M(1,M3 
CONTINUE 
GO  TO  1G1 
MBG=MRG+1 
MBGS(MBG)=MC 
CONTINUE 
IF  (MBG.EQ.O 
MA=MA+l 
MAMBIG(MA,J  ) 
MAMBIG(MA,2 ) 
M£MBIG(MA,2) 
MAMBIG(MA,4) 
MAMBIG(MAt5) 
DO  104  M35=l 
MAMBIG(MA, p+ 
KA=MBG+5 
GOTO  90 
MU=MU*1 
MUNI0(MU,1 ) = 
MUNI0(MU,2)= 
MUNIO( VU,3)= 
MUNI0(WU,4)= 
CONTINUE 
PINO  UNIQUE  CONN 
ENDNODES  OF  A  P 
NSIMP=0 
NSEP=0 

DO  IC7  M39=l 
MCP=C 

M40=TRUNK(M8 
DO  1C8  M41=l 
CHECK  =  C. 
RELP0S=1 
MEISER=0 
MCOL=MAINBP( 

MUFLA=C 
MCOT=C 
MC0T1=0 
M42=MSER(M41 
DO  1C9  M43  =  l 
IF(MSLET  (M4 
M44=MSLET(M4 
IF  (M44.GT.1 
MC0T1=MC0T!+ 
DC  11C  M45  =  l 
IF  (CHECK. EO 
CONTINUE 
DO  228  N75  =  l 
N76=TRUNK(M3 

IF  (N76.EQ.M 

CONTINUE 

IF(MBRT(M41, 

CHECK*! 

IF  (N75.LT. M 

GOTO  229 

CONTINUE 

CONTINUE 

IF(MBRT(M41 , 

MCOT=MCOT-»-l 

IF(MSLFB(MA1 


M26-1) )G0  TO  100 


)=MBR8(M10f M23,M25) 


0.0P.MC(2).EQ.0)G0  TO  1  >5 


tM31 
,  M32 
3>.EQ.MCOM(2,M34) )  GO  TO  103 


0M(1 ,M33) 

>  GOTO  105 

=F(1 , MCOL) 

=  MRG 

=  NEX 

=  M10 

=  1 

,MRG 

M35 >=MBGS(M35) 


F(l ,MCOL) 

NEX 

Mir 

1 

ECTION  TO  ONE  OF  THE 
ATH  OF  TRUNK  BRANCHES 


tTR 

,M?9) 
,MBR 


M41) 


J 

,M42 

lt^43).E0.r | GO  TO  109 

1 ,M43) 

)      MEISER=1 

I 

,  M44 

.1)  GOTO  229 


2?8 
229 


tTC 

,N75) 

40)  GOTO  228 

M43tM45l.NE.N76)  GOTO  228 

39)  RELP0S=2 

M43fM45).NE.M40 )G0  TO  110 
,M43I.EQ.C)  MUFLA=1 
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110  CONTINUE 
109   CONTINUE 

IF  (MEISER.EQ.O)  GOTO  108 
CONTINUE 

IFfMCOTl.EQ.MCOT.AND.MCOT.NE.fMGO    TO    111 
GO    TO    7.C8 

111  MCP  =  MCP«-1 

IF    (MUFLA.EO.l)    GOTO    112 

MC0MP(MCPI=M41 

GO    TO    108 

112  CONTINUE 
MU=MU+1 
MCP=0 

MUMI0(MU,1)=F(1 ,MCOL) 
MUNI0<MU,2)=M40 
MUNI0(MU,3)=M41 
MUNIQ(MU,4)=RELP0S 

108   CONTINUE 

IF  (MCP.EQ.O)  GOTO  107 

CONTINUE 

DO  113  M47=1,MCP 

MC1=0 

MBG=0 

M48=MC0MP(M47) 

MA9=MSER(MAB) 

IF(M49.NE.1)G0  TO  114 

MC1=MSLEB(M48,1 ) 

DO    115    M36=1,MC1 

115  MC0M(1.M36)=MBPR(M48,! ,M36) 
GO    TO    122 

114      CONTINUE 

DO    116    M5C=ltM49 

IF    (MSLET(M48,M50).E0.0)    GOTO    116 

M51*MSLET(M48,M50) 

DO  117  M52=lfM51 

IF(MBRTf M48,M49,M52».EO.M40)GO  TO  118 

117  CONTINUE 
GO  TO  116 

118  MBG=MBG«-1 
MBGS(MBG)=M50 

116  CONTINUE 
M54=MBGS(1 ) 
M55=MSLEB(M48,M54) 
DO  119  M56=l,M55 
MCOT=C 

M57=MBRB(M48,M54,M56) 
DO    120    M58=2,MRG 
M59=MBGS(M53) 
M6C=MSLEB<M48,M59) 

DO  121  M61=1,M6C 

IF(MBRB(M46,M59, M61 I.NE.M57)  GO  TO  121 

MCOT-MCOT+1 

121  CONTINUE 
120      CONTINUE 

IF(MCOT.NE. (MBG-1 ) ) GO    TO    119 

MC1  =  MCH-1 

MC0M(1,MC1)=M57 

119  CONTINUE 
IF(MC1.NE.C)G0  TO  122 
MU=MU+1 

MUNI0(MU,1)=F(1 ,MCOL) 
MUNI0(MUt2)=M4C 
MUNI0(WUf3)=M41 
MUNI0(MU,4)=3 

IF  (M39.NE.1.AND.M39.NE.TR)  MUN I0( MU, 4) =R ELPOS 

CONTINUE 

GOTO  113 

122  MA=MA*1 
MC0L=MAINBP(M48) 
MAMBIG(MAfl )=  F(1,MC0D 
MAMBIG(MA,2)=MC?. 
MAMBIG(MA,3)=M40 
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MAMBIG(MA, 4)=M48 

MAMBIG(MA,5)=3 

IF  (M39.NE.1.AND. 

CONTINUE 

DO  123  M62=1»MC1 

123  MAMBIG(MA,f  4-M62)  = 
KA=MCl+5 

113   CONTINUE 
107   CONTINUE 
MGR=0 

124  CONTINUE 

C  STORE  CONNECTIONS  IN 
IF  (MU.EQ.M  GOTO 
CONTINUE 

IF  (MA.EQ.C)  GCTO 
MCT=0 

DO  125  M65=1,MU 
IF(MCT.EO.MA)GCTO 
CONTINUE 
DO  127  M66=1,MA 
IF(MAMBIG(M66, i). 
M67=MAMBIG(M66,2) 
DO  128  M68=6,M67 
IF  (MAMBIG( M66,M6 

128  CONTINUE 
GO  TO  127 

130  MAMBIG(M66,l)=r 
MCT=MCT+1 

127  CONTINUE 
12  5  CONTINUE 
126   CONTINUE 

DO  131  M69=1,MU 

MGP=MGR+1 

GRAPH(MGR,1 )=MUNI 

GPAPH(MGR,?)=MUNI 

GRAPH(MGR.3)=MIJN! 

M70=MUNIQ(M69,3) 

ME=ME+1 

MEND(ME)=MUNIQ(M6 

131  MSEP(M7C)=-MSER(M 
133   MU=0 

IF  (MA.EQ.M  GOTO 
CONTINUE 
DO  134  M71=1,MA 
IF(MAMBIG(M71,1). 
CONTINUE 

DO    141    MAC=1,MGR 
IF     (GRAPH<MAC,1). 
141       CONTINUE 
MGR=MGR+1 
GRAPH(MGR,1)=MAMB 
GRAPH(MGR,2)  =  W>AMB 
GRADH(MGR,3)=MAMB 
M72=MAMBIG(M71,4) 
ME=ME+1 

MEND(ME)=MAMBIG(M 
MSER(M72>=-MSER(M 
M73=MAMBIG(M71,2) 
DO  136  M74=6,M73 
DC  129  MAK=1,MGR 
IF  (GRAPH(MAK, 1). 

129  CONTINUE 
MGP=MGR+1 
GRAPH(MGR,1 )=MAMB 
GPAPH(MGR,2 )=GPAP 
GRAPH  <MGR,3)=7 
DO  137  M75=1,MBR 
M76=MAINBR<M75) 
IF  <F(1,M76).NE.M 
MSER(M75)=-MSEP(M 
ME=ME+1 
MEND(ME)=MAMBIG(M 


M39.NE.TR)  MAMBIG( MA ,5 ) =RELPOS 


MC0M(1,M62) 


GRAPH 
133 

126 


126 


EO.C )G0  TO  127 
+  5 

8).FQ.MUNIQ(M65,m  GOTO  130 


3(M69,1) 

0(M69,2) 
Q(M69,4) 


9Xl) 
70) 


138 


EQ.OGO  TO  134 


E0.MAMBIG(M71t! ))  GOTO  134 


IG(M71f 1) 

IG(M71,3) 
IG(M71,5) 


71,1) 
72) 
+  5 


FQ.MAM8IG(M71,M74) )  GOTO  136 


IG(M71,M74) 
H(MGR-1,1) 


AMBIG(M71,M74) )GOTO  137 
75) 

71,H74) 
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137  CONTINUE 
136  CONTINUE 
134   CONTINUE 

138  MA=0 

139  NSIMP=^ 

140  CONTINUE 

C  FIND  CONNECTIONS  TO  ALREADY  CONNECTED 
C  MAINBRANCHES. IF  CONNECTION  IS  UNIQUE 
C   NSIMP  IS  ZERO 

CCT=0 

DO  142  N1=1.ME 

N3=MEND(ME+1-N1 > 

MCP=0 

DO  143  N4=1 ,MBP 

IF  (MSFR(NM.LE.O)  GOTO  143 

MCOTl=0 

MCOT=0 

N5=MSER(N4) 

DO  144  N6=1,N5 

IF  (MSLEB(N4,N6>.E0.0)  GOTO  144 

MC0T1  =  MC0TH-1 

N7=MSLEB(N^,N6) 

DO  145  N8=l ,N7 

IF  (MBRB(N4,N6,N8).NE.N3)  G0T0145 

MC0T=MC0T*1 

IF  (NSIMP. FQ.1C)  GOTO  146 

CONTINUE 

IF  (N7.EQ.1)  GOTO  147 

145  CONTINUE 
144   CONTINUF 

IF  (MCOTl.FQ.MCOT.AND.MCOT.NE.'M  GOTO  146 

CONTINUE 

IF  (MCOT.NF.C)  GOTO  146 

GOTO  143 

146  MCP=MCP*1 
MCOMP(MCP)=N*> 
GOTO  143 

147  CONTINUE 

IF  (CCT.EO.O)  GOTO  2J0 
DO  209  C1=1,CCT 

IF  (N4.E0.C0NECT(C1))  GOTO  211 
239   CONTINUE 

210  CONTINUE 
CCT=  CCT  +  1 
C0NECT(CCT|=N4 
MU=MU+1 
N8=MAINBR(N4) 
MUNIQ(MU, 1)=F( 1,N8) 
MUNI0(MU»2)=N3 
MUNI0(MU,3)=N4 
MUNIQ(MU,4)=3^NSIMP 

211  CONTINUE 
143   CONTINUE 

IF(MCP.EQ.C)  GOTO  142 

CONTINUE 

DO  148  N9=1,MCP 

MC1=0 

MBG=0 

N10=MC0MP(N9) 

MC0L=MAINBR(N10) 

N11=MSER(N10) 

IF    (Nll.NE.l)    GOTO    150 

MC1=MSLEBIN1C,1 ) 

NIC  =  0 

DO    149    N12=i,MCl 

IF(MBRB{N1C,1,N12).EQ.N3)    GOTO    149 

NIC=NIC*1 

MC0M(lfNIC)=MBRB(N10,l,N12) 

149  CONTINUE 
MC1=NIC 
GOTO    159 

150  CONTINUE 
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00  151  N13=ltNll 

IF  (MSLEB(N10,N13).EQ.O)  GOTO  151 
N14  =  MSLEB(NKtN}.3) 
DO  152  N15=1,N14 

IF  (MBRB(N10, N13,N15) .EQ.N3)  GOTO  153 
15?   CONTINUE 
GOTO  151 

153  MBG=MBG*1 
MBGS(MBG)=N13 

151   CONTINUE 

IF  (MBG.EQ.l)  GOTO  161 

N16=MBGS<1) 
N17=MSLEB(NU ,N16) 

DO  154  N18=1,N17 

MCOT=0 

N19=MBRB(N1^,N16,NI8) 

IF  (Ni9.E0.N3)  GOTO  154 

CONTINUE 

DO  155  N20=1,MRG 

N21=MBGS<N20) 

N22=MSLEB(N10,N21) 

DO  156  N23=l.  ,N2? 

IF(MBRB(N19,N22 , N23 ) . NE.N19 )  GOTO  156 

MC0T=MC0T+1 
155   CONTINUE 
155   CONTINUE 

IF(MCOT.NE. (MBG-1 ) )  GOTO  154 

MC1=MCT+1 

MC0M<lfMCll=N19 

154  CONTINUE 
IF(MCl.EQ.C)  GOTO  161 

159  CONTINUE 

IF  (CCT.EO.O)  GOTO  212 
DO  213  C2=1,CCT 
IF(N10.E0.C0NECT(C2))  GOTO  148 

213  CONTINUE 
212   CONTINUE 

CCT=    CCT+1 

C0NECT(CCT)=N1C 

MA=MA+1 

MAMBIG(MA,1 )=FQ,MCOL) 

MAMBIG(MA,2)=MC1 

MAMBIG(MAt3)=N3 

MAMBIG(MA,4)=N1C 

MAMBIG(MAf5)=3+NSIMP 

DO  160  N24=1,MC1 

160  MAMBIG(MA,5*N24)=MC0M(lfN24) 
KA=MCl+5 

GOTO  148 

161  CONTINUE 
MU=MU«-l 

IF  (CCT.EO.O)  GOTO  214 
CONTINUE 
DP  215  C3=1,CCT 

IF(N10.E0.C0NECT(C3))  GOTO  148 
215   CONTINUE 

214  CONTINUE 
MUNI0(MU,1)=F(1 ,MCOL) 
MUNIQ(MU,2)=N3 
MUNI0(MUf3)=NlC 
MUNI0(MU,4)=3+NSIMP 

148  CONTINUE 
142   CONTINUE 

IF(<mJ+MA).NE.C)  GOTO  124 

MCOUNT=0 

DO  163  N25=1,MBR 

IF(*SER(N25).LE.G)  GOTO  173 

MC0L=MAINBR(N25) 

IF(NSIMP.NE.O)  GOTO  164 

NSIMP=10 
C  NO  UNIQUE  CONNECTION  COULD  BE  MADE 

GOTO  140 
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C  CONNECT  SO  AS  TO  FORM  A  LINEAR  PATH  WITH 
C   MAIN  BRANCHES 
164   N29=MSEP(N25) 

DO  166  N26=1,MGR 

N27=MGR*1-N26 

N28=GRAPH(N27,1) 

00  168  N30=ltN29 

IF<MSLEB<N25,N3C).EQ.C)GOTO    168 

N31=MSLEB(N25,N30) 

DO  169  N32=lfN31 

IF(MBRB(N25,N30,N31).NE.N28)  GOTO  169 

MU=MU«-1 

MUNI0(MU,1)=F(1 ,MCOL) 

MUNlO(MUf2)=N28 

MUNI0(MU,3)=N25 

MUNI0(MU,4)=3*NSIMP 

GOTO  126 

169  CONTINUE 
168  CONTINUE 
166      CONTINUE 

C    CONNECT    SO    AS    TO    FORM    A    LINEAR    PATH    WITH 
C       TRUNK    BRANCHES 

DO    170    N34=1,TR 

N38  =  TRUNMM3,N34) 

DO  171  N35=i,N29 

IF(MSLET(N25,N35).EQ,C)GOTO  171 

N36=MSLET(N25,N35) 

00  172  N37=1,N36 

IF(MBRT(N2,:,N35,N37)#NE.N38)  GOTO  jl  72 

MU=MU*1 

MUNIO  (MU,1 )=F(1 ,MCOL) 

MUNIO  (MU,2)=N38 

MUNIO  (MU,3)=N2  5 

MUNIQ(MU,4)=H-NSIMp 

GOTO  126 

172  CONTINUE 
171   CONTINUE 

170  CONTINUE 

C  MAIN  BRANCH  STARTS  A  SEPARATE  PART 
MU=MU+1 

MUNI0(MU,1)=F(1 ,MCOL) 
MUNI0(MU»2)=0 
MUNIO(MU,3)=N25 
MUNI0(MU,4)=6 
GOTO  126 

173  MC0UNT=MC0UNT*1 
IF(MCOUNT.EO.MBR)    GOTO    216 

163  CONTINUE 
216   CONTINUE 

IF  (L.EO.C)  GOTO  175 

CONTINUE 

CALL  LIMBS 
175   CONTINUE 

CALL  LINKS 

CONTINUE 

CALL  OUTPUT 
89    CONTINUE 

350  FORMAT( »1 • ,T30, 'F-MATRIX'/) 

351  F0RMAT(//T30, • F-MATRI X, PART ITIONED' / ) 

353  F0RMAT(//1H  , T30 , « MAINLOOP  IS  ROW  *  »,I2) 

354  FORMAT  ( //26X, 12 , T3C , 'TRUNKBRANCHES • ) 
309   FORMATUH  .27X.2CI3) 

355  FORMAT   ( /26X, 12 , T30 ♦ 'MAI NBRANCHES •> 

356  FORMAT   ( /26X, 12 , T30 , » L IMBS • ) 

357  FORMAT(//T30f 'COMMON  NODE  REQUIREMENTS') 
353   FORMAT   ( //26X, 12 ,T3C ,• TRUNKS • ) 

STOP 
END 

SUBROUTINE  OUTPUT 
C  THIS  SUBROUTINE  CONVERTS  NUMBERCODE  FOR 
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CONNECTIONS  INTO  WORDS 

IMPLICIT  INTEGER  (A-Z) 

01  MENS  I  ON    F(30f3<M  tLIMB(2CI  ,NOTLIM(  2C  )  ,MAINBR(20) 

DIMENSION    truNRP(2C) 

DIMENSION    MEND(  30)fGRAPH(6C  ♦  ?>  t  TRUNK  (20,2  C  ),Ft>ART(3C) 

COMMON/ LLO/    ROW  ,  COL  ,  F  ,TR,  TRUNRP ,  M8  ,  TRtJNK  ,  MBR  ,  MA  IN  BR  , 
1MEND,MGR,GPAPH,L,LIMB,ME 

WRITE(6,359)     MB 

WRITE(6,3C9)     (TRUNK(M8,  I)  ,I  =  1,TP  ) 

WRITE(6,36D 

IF(MBR.NE.C)     WRITE(6,343) 

M0THER=10 

DO    314    N63=T,MGR 

MOTHER-MOTHER*! 

IF(N63.EQ. (MBR+1)  )     WRITE(6,34M 

CONTINUE 

IF(N63.E0.  MBR  +  L  +  1)  )WRITE(6,345) 

CONTINUE 

IF(N63.E0.  (MBR+L+1)  )M0THER=-1 

CONTINUE 

IF(  MOTHER,  r-Q.O)    GO    TO    341 

CONTINUE 

IF( MOTHER. GT.1C)  GOTO  361 

MOTHER*- 1. 
361   WPITE(6,31C)GRAPH(N63,1 > 

GOTO  346 
341   WRITE(6.342) 
346   LOKT=GRAPH(N63,3) 

LOCT=LOKT 

IF  (L0KT.GT.7)  LCCT=LOCT-iC 

CONTINUE 

IF  (LOCT.EO.l)  GOTO  321 

CONTINUE 

IF  (L0CT.E0.2)  GOTO  322 

CONTINUE 

IF  (L0CT.EQ.3)  GOTO  323 

CONTINUE 

IF  (L0CT.EQ.4)  GOTO  324 

CONTINUE 

IF  (L0CT.E0.5)  GOTO  325 

CONTINUE 

IF  (L0CT.EQ.6)  GOTO  326 

CONTINUE 

IF  (L0CT.EQ.7)  GOTO  327 

321  WRITE(6,33l) 
GOTO  316 

322  WPITE(6,332> 
GOTO  316 

323  WRITE(6,333) 
GO  TO  316 

324  WRITE(6,33M 
GO  TO  316 

325  WPITE(6,335) 
GOTO  316 

326  WRITE(6,336) 
GOTO  314 

327  WRITE(6,337) 
WRITE(6,3ie)  GRAPH(N63-1,1) 
WPITE(6,33G) 

GOTO  314 

316  CONTINUE 
WRITE(6,317)  GRAPH(N63,2) 

IF  (L0KT.GT.7)  WRITE  (6,338) 

314  CONTINUE 

315  FORMAT    (/31X,I2) 

317  F0RMAT(1H*,52X,I2) 

318  FCRMAT(lH-»-,49X,I2) 

331  FORMAT( •♦• ,T39, 'LEFT  OF') 

332  FORMAT( •  +  • ,T39, 'RIGHT  OF  •  > 

333  FORMAT( '♦• ,T39, 'OUTSIDE  OF') 

334  FORMAT(  •♦•  rT39,  •  INSIDE  OEM 

335  FORMAT( •♦» ,T39, 'EITHER  SIDE  OF') 
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309  FORMATdH  t27X,20I3) 

336  FORMAT  (  •♦•  tT39,  'STARTS  A  SEPARATE  PART') 

337  FORMAT! •♦• ,T39, 'AND  BRANCH*') 

339  FORMAT ( •♦• ,T54, • INTERCHANGEABLE  •  ) 

338  FORMAT( '♦• ,T61 ♦ • (NOT  UNIQUE)1! 

342  FORMAT(T3,'  •) 

343  FORMAT(/T3C,'MAIN: • ) 
300  FORMAT  (1H  ,28X,30I2) 

344  FORMAT(/T3ff 'LIMB: • ) 

34  5  FORMAT(//T30,'LINK: • ,T39 ,  • FROM/ tqi t T5 3, • BRANCH' / ) 

359  FORMAT(//T30, 'GRAPH  FOR  TRUNK • , I  2 ♦ J H:  ) 

36C  FORMAT(//T30, 'BRANCH' ,T39, 'SIDE  OF • ,T53 ,' BRANCH' ) 

RETURN 

END 

SUBROUTINE  LINKS 
C  THE  POSITION  RELATIVE  TO  THE  TREE  IS  FOUND. 
C   THE  CONNECTION  PATTERM  IS  INDICATED  BY  CASF 
IMPLICIT  INTEGER  (A-Z) 

DIMENSION  F(30,30> , LI  MB (21 ) ,NOTLIM(2C  ),MAINBR(2C  ) 

DIMENSION  TRUNBR(20) 

DIMENSION  MEND(30) , GRAPHt 60 , 3 ) , TPUNK ( 20 ,2C  ),FPART(3>  I 

COMMON/ LLO/  ROW , COL , F , TR , TRUNBR , M8 , TRUNK , MBR , MA  I NBP , 
1MEND.MGR,GRAPH,L,LIMB,ME 

TGR  =  C 

DO  183  N5C=2,ROW 

TSUM=0 

IF  (TR.EQ.C)  GOTO  184 

CONTINUE 

DO  185  N51=1.TP 

N52  =  TRUNBR(N5.\  ) 

IF  (F(N50,N52).EQ.l)  TSUM=TSUM+1 

185  CONTINUE 
134   MSUM=0 

IF  (MBR. EG. 0)  GOTO  186 

DO  187  N53=1,MBR 
N54=MAIN8R(N53> 

IF  (F(N50,M  541.E0.1)  MSUM=MSUM+i 
137   CONTINUE 

186  LSUM=0 

IF  (L.EO.CI  GOTO  188 

DO  189  N55=1.L 

N56=LIMB(N55) 

IF  (F(N50,N56).E0.1>  LSUM=LSUM*1 
189   CONTINUE 
188   CASE=6 

IF  (LSUM.NF.O)  CASE=CASE-3 

CONTINUE 

IF  (MSUM.EO.O.AND.TSUM.EO.O)  CASE=CASE-3 

CONTINUE 

IF  (MSUM.NF.O)  CASF=CASE-1 

CONTINUE 

IF    (TSUM.NE.O.AND.MSUM.NE.O)    CASE=CASE-1 

CONTINUE 

IF  (CASE.E0.6)  GOTO  196 

CONTINUE 

IF  (CASE.NE.4)  GOTO  230 

NT=0 

DO  227  NS1=1.MBR 

NS2=MAINBR(NS1) 

IF  (F(N50,  NS2).EQ.O)  GOTO  227 

CONTINUE 

DO  224  NS3=1,MGR 

IF  (GRAPH(NS3,2).EO.F(l,NS2) )  GOTO  227 

224  CONTINUE 
NT=NT*1 

IF  (NT.EO.2l  GOTO  225 

NTMB=F(1,NS2) 

GOTO  227 

225  MGR=MGR*1 
GRAPH(MGR,1»=F(N50,1) 
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GRAPH(MGRt2)*NTMB 

GRAPH(MGR,3  )=3 

MGR=MGR*1 

GRAPH(MGR,1 )=F(N50.1) 

GRAPH(MGR,?)=F(1,NS2) 

GRAPH(MGR,3>=3 

GOTO  183 
227  CONTINUE 
230   CONTINUE 

IF  (CASE.GE.4)  GOTO  190 

CONTINUE 

DO  191  N57=1,L 

N58=LIMB(L+1-N57) 

IF  (F(N5O,N58).EQ.0>  GOTO  191 

DO  192  N70=1,MGR 

IF  (GRAPH(MGR+1-N7G,1).NE.F(1 ,N58) )  GOTO  192 

MG=MGR+1 

GRAPH(MG  ,1I=F(N50,1> 

GRAPH(MG  ,2)=M1,N58) 

GRAPH(MG  ,3)=3 

GOTO  22  6 

192  CONTINUE 
191  CONTINUE 
226   MGR  =  MGR«-1 

IF  (CASE.EO.l)  GOTO  194 
CONTINUE 

IF  (CASE. EG. 2)  GOTO  190 
CONTINUE 
DO  193  N71=1,L 
N72  =  LIMB(N"'\) 

IF  (F(N5C,N72>.EQ.0)  GOTO  193 
CONTINUE 

DO  195  N69=1,MGR 

IF  (GRAPH(MGP+1-N69,1 I.NE.FIl ,N72) )  GOTO  195 
CONTINUE 

IF  (GRAPH(MGR+1-N69,3).EQ.7.AND.CASE.EQ.3>  GOTO  195 
CONTINUE 

IF  (CASE.E0.3»G0T0  197 
MG=MGR+1 

GRAPH(MG  ,1 )=F(N50,1) 
GRAPH(MG  »2)=F(1,N72) 
GRAPH(MG  ,3)=3 
GOTO  197 
195   CONTINUE 

193  CONTINUE 

197  TGR=GRAPH(MGR+1-N69,2) 
IF     (CASE.EQ.O)     MGR  =  MGRU 
CONTINUE 

IF    (CASE. EG. 0)     GOTO    183 
IF    (GRAPH(MGR*l-N69,3).EQ.l)    P0S0NE=2 
CONTINUE 

IF     (GRAPH(MGR+1-N69,3».E0.2)     P0S0NE=1 
GOTO    196 
190       MLOW=100 
MHIGH=0 

DO    198    N59=1,MBP 
N6C=MAINBR(N59) 
IF    (F(N50,N60).EQ.G)    GOTO    193 
CONTINliE 
00    199    N61=l .ME 

IF     (MEND(N61).NE.F(1,N60) )     GOTO    199 
CONTINUE 

IF     (N61.LT.M10W)     ML0W=N6l 
CONTINUE 
IF     (N61.GT.MHIGH)     MHIGH=N61 

199  CONTINUE 

198  CONTINUE 

IF  (CASE-4)  200,201,2CC 

200  MGR=MGR+1 
GRAPH(MGR,1 )=F(N50,1) 
GRAPH(MGR,3»=4 
GRAPH(MGRf2)=MEND(ML0Wl 
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IF  (CASE.E0.2)  GOTO  183 

201  MGR=MGR+1 
GRAPH(MGR,1 »=F(N50,l) 
GRAPH(MGR,3)=3 
GRAPH(MGR,2)=MEND<MHIGH) 
IF  (CASE.E0.5)  GOTO  183 

194   CONTINUE 

DO  202  N62=1,MBR 

N63=MAINBR(N62) 

IF  <F(N50,N63).EQ.O)  GOT02C2 

DO  203  N64=1,MGR 

IF  (F(l,N63).NE.GRAPH<N64,l ) )  GOTO  203 

CONTINUE 

DO  204  N65=1,TR 

IF(GRAPH(N64,2) .NE.TRUNKCM8 »N65) )  GOTO  204 

TGR=TRUNK<M8,N65) 

IF(GRAPH(N64,3).EQ.1>  P0S0NE=2 

CONTINUE 

IF(GRAPH(N64,3).EQ.2)  P0S0NE=1 

GOTO  196 

204  CONTINUE 
203   CONTINUE 

202  CONTINUE 
196   MLEFT=50 

IF  (CASE.EQ.6)  TGR=0 

MRIGHT=0 

TPOS=0 

DO  205  N66=i,TR 

N67=TRUNBR(N66) 

IF  (F(N50,N67).EQ.O)  GOTO  205 

CONTINUE 

DO  206  N68=1,TR 

IF  (TRUNK(M3,N68).EQ.TGR)  TP0S=N68 

CONTINUE 

IF  (TRUNK(M8,N68) .NE.F(1,N67) )  GOTO  2T6 

CONTINUE 

IF  (N68.LT. MLEFT)  MLEFT=N68 

CONTINUE 

IF  (N68.GT.MRIGHT)  MRIGHT=N63 

206  CONTINUE 

205  CONTINUE 

IF  (CASE.E0.6)  GOTO  2C7 

CONTINUE 

IF  (TPOS.EO.MLEFT)  GOTO  208 

207  MGR=MGR+1 
GRAPH(MGR,1 )=F(N50tl) 
GRAPH(MGR,3)=1 

IF  (MRIGHT.EO.MLEFT.AND.MLEFT.EO.TPOS ) 
1  GRAPH(MGR,3)=P0S0NE 
GRAPH(MGR,2)=TPUNK(M8,MLEFT) 
IF  (CASE.NE.6)  GOTO  183 

208  MGR=MGR-H 
GRAPH(MGR,1 )=F(N50,1) 
GRAPH(MGR,3)=2 

IF  (MRIGHT.EQ. MLEFT. AND. MLEFT. FO. TPOS ) 
1  GRAPH(MGR,3)=POS0NE 
GRA  PH ( MGR , 2  I =T RUNK ( M8 , MR  I GHT ) 
133   CONTINUE 
RETURN 
END 

SUBROUTINE  LIMBS 
C  THE  LIMBS  ARE  CONNECTED  SO  AS  TO  FORM  LINEAR 
C   PATHS  WITH  MAIN  OR  TRUNK  BRANCHES 
IMPLICIT  INTEGER  (A-ZI 

DIMENSION  F(30,30),LIMB(2C) ,NOTLIM ( 20 ) ,MA INBR( 2C) 

DIMENSION  TRUNBRC20) 

D I  MENS  ION  MEND( 30 ) , GRAPH( 60 , 3 ) , TRUNK ( 20 , 2C ) . FPART ( 30 > 

COMMON/ LLO/  ROW, COL ,F , TR, TRUNBR,M8 , TRUNK , MBR , MA INBP , 
1MEND,MGR,GRAPH,L,LIMB,ME 

CONTINUE 
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DO  177  N42=2fROW 

ROWLIM=0 

00  176  N40=1»L 

N43=LIMB(NV;  ) 

IF  <F(N42,N4n.FQ.0l  GOTO  3  76 

P0WLIM=R0WLIM+1 

IF  (ROWLIM.EQ.I JGOTO  219 

MGR=MGR*1 

GRAPHCMGR,! )  =F  ( 1 ,N41) 

GRAPH(MGR,2 ) =GRAPH ( MGR-1 t 1 ) 

GRAPHCMGR f? )=7 

GOTO  176 
219   NOMBR=rt 

MSTORS=C 

DO  178  N43=1,MBP 

N44*MA1NBR(N43I 

IF  <F(N42,N44).EO.O)  GOTO  178 

N0MBR=1 

DO  179  N45=1,ME 

IF  (F(l,N4M.NE.MEND<N45n  GOT03  79 

CONTINUE 

IF  (N45.LT. MSTORS)  GOTO  179 

MST0RS=N45 
179   CONTINUE 
178   CONTINUE 

IF  (NOMBR.EO.O)  GOTO  18C 

MGR=MGR+1 

GRAPH (MGR,1 >=F(] ,  N41) 

GRA PH ( MGR» 2 »=MEND( MSTORS) 

GRAPH(MGR,?)=3 

GOTO  176 
13C   MST0RT=50 

NOTTBR=0 

DO    181    N46=1,TR 

N47=TRUNBR(N46) 

IF  <F(N42,N47>.E0.0)  GOTO  181 

CONTINUE 

DO  182  N48=1,TR 
IF  <F<1,N^7).NE.TRUNMM8,N43))  GOTO  182 

N0TTBR=1 

IF  (N48.GT.MST0RT)  GOTO  182 

MST0RT=N48 
132   CONTINUE 
191   CONTINUE 

MGR=MGR+1 

GRAPH(MGR,1»=F(1,N41) 

GRAPH(MGR,?l=0 

GRAPH(MGR,3)=6 

IF  (NOTTBR.EO.C)  GOTO  176 

GRAPH(MGR,2)=TRUNK(M8,MSTORT) 

GRAPH(MGR, 31=11 

176  CONTINUE 

177  CONTINUE 
RETURN 
FND 

SUBROUTINE  CHAINS(K) 
C  CHAINS  ARE  CHECKED  FOR  DUPLICATIONS 
IMPLICIT  INTEGER  (A-Z) 
DIMENSION  KSTOPE  ( 20 ) .KCHAI N< 23 , 10 , 20 ) , KC ( 10 ) 
COMMON  /CH/  KCHAIN»KSTOREtKC 
IF  (KC(K).EO.O)  GOTO  1 
KR=KC(K» 
DO  2  K1=1,KR 
KSAME=C 
DO  3  K2=1.K 
K3=KCHAIN(K,K1,K2I 
DO  3  K4=1,K 

IF  (K3.E0.KSTORE(K4) )  KSAME=KSAME*I 
3     CONTINUE 

IF(KSAME.EO.K)  GOTO  5 
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2  CONTINUE 

1  KC(K)=KC(K>«-1 
KR=KC(K) 

DO   4   K3=1,K 

4  KCHAIN  (K,KR,K3)=KSTORE(K3) 

5  RETURN 

END 

SUBROUTINE  JTRONE  (N,T,M 
C  DELETE  THOSE  JTRCH'S  WHICH  HAVE  BEEN  COMBINED, 
C   ALSO  DELETE  THEIR  VARIATIONS 
IMPLICIT  INTEGER  (A-Z) 

DIMENSION  JTRCH(30,30),LENGTH(20)  ,JLAST(10),JFIRST(13) 

DIMENSION  JTFMP  (30) 

DIMENSION  JADD(30,9) 

COMMON/PT/JTPCH,LENGTH,JTWO,JLAST,JFIRST, J ADD, 
1JTEMP  ,JT,JF.JL 

IF  (N.EO.O)  GOTO  1 

DO  2  J1=1,N 

IF  (IABS(LENGTH( Jl) ).NE. I ABS ( LENGTH ( T ) ) )  GOTO  2 

J2=IABS(LENGTH(T) ) 

DO  3  J3=1.J2 

IF  (JTRCHlJl ,J3).NE«JTRCH(T, J3) )  GOTO  1 

3  CONTINUE 

2  CONTINUE 
1     RETURN 

4  RETURN1 
END 

SUBROUTINE  PERMUT  (J, NEW,*) 
C  GENERATE  VARIATIONS  OF  NEW  JTRCH'S  SUCH 
C   THAT  PERMUTABLE  BRANCHES  APPEAR  ONCE  IN 
C    FIRST  OR  LAST  POSITION  OF  A  JTRCH 
IMPLICIT  INTEGER  (A-Z) 

DIMENSION  JTRCH (30, 30) , LENGTH* 20) ,JLAST(10),JFIRST(r) 

DIMENSION  JTEMP  (30) 

DIMENSION  JADD(30,9) 

COMMON/ PT/ JTRCH, LENGTH, JT WO, JL AST, J  FIRST, J ADD, 
1JTEMP  ,JT,JF,JL 

IF  (J.EO.l)   GOTO  12 

CONTINUE 

IF  (J.E0.2)   GOTO  11 

CONTINUE 

IF  (J.FQ.3)   GOTO  3 

CONTINUE 

IF  (J.E0.4)   GOTO  10 

CONTINUE 

IF  (J.EQ.5)   GOTO  5 

CONTINUE 

IF  (J.EQ.6)   GOTO  6 

3  JTURN=JFIRST(1) 
JF=0 

JL=1 

JLAST(1)=JTURN 
GOTO  8 

5  JTURN=JLAST(2) 
JF  =  1 

JL  =  1 

JFIRST(1)=JTURN 
GOTO  8 

6  JTURN=JFIRST(2) 
JF  =  1 

JL  =  1 

JLAST(1)=JTURN 
3     J1=LENGTH( JTURN) 

IF  (Jl.LT.C)  GOTO  10 

J2=Jl/2 

DO  9  J3=1.J2 

JTP=JTRCH( JTURN, J3) 

JTRCH ( JTURN, J3)= JTRCH (J TURN, J 1-J3  +  1  ) 
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9     JTRCH( JTURN,J1-J3*1)=JTP 
U    CONTINUE 

IF  (JT.NF.C)  GOTO  11 

NEW=C 

GOTO  16 

11  CONTINUE 

IF  (JT.GT.ll  GOTO  12 
NEW  =  1 

JADD(1,1)=JTEMF(1) 
GOTO  16 

12  NEW=0 

DO  13  J4=1,JT 

IF  (J   .EO.l.AND. ( ( J4+1 )/2).EO. (J4/2) )  GOTO  13 

NEW=NEW*1 

JCOUNT=0 

DO  14  J5=1,JT 

JR=J4+J5-1 

IF  ( JR.GT. JT) JR=JR-JT 

JCOUNT=JCOUNT+l 

JADD(NEW,JCOUNT)=JTEMPUR) 

IF    ( J.EQ.1.0R. JTWO.EQ.O)    GOTO    14 

JNEW=NEW*-JT 

JADD( JNEW,JCOUNT)=JTEMP(JT>l-JR> 

14  CONTINUE 

13  CONTINUE 

IF  (J.EQ.l)  GOTO  15 

CONTINUE 

IF  (JTWO.EO.C)  GOTO  16 

NEW=JNEW 

GOTO  IS 

15  RETURN 

16  RETURN! 
ENO 

SUBROUTINE  NODUP 
C  DELETE  THOSE  OF  JTPf.H'S  WHICH  ARE  DUPLICATED 
C   OR  WHICH  VIOLATE  A  COMMON  NODE  REQUIREMENT 
IMPLICIT  INTEGER  <A-Z) 

DIMENSION  JTRCH(30,3O) , LENGTH (2^) , J  LAST ( 10 ) , JF  IRS T( 20) 

DIMENSION  JTPMF  (30) 

DIMENSION  MBOTH(20)  ,MNO(20l  ,MLORF(2")) 

DIMENSION  JADD(30,9) 

COMMON  /NOD/JDELET (40l,JD,COMB<4O),C, J 

COMMON/PT/JTRCHtLENGTHt JTWOtJLASTf JFIRST, JADO, 
1JTEMP  ,JT,JF,JL 

C=0 

MB=C 

MO=0 

MLF=C 

DO  1  M1=1,J 

GO=0 

MLE=LENGTH(M1) 

DO  2  M2=1,JT 

IF( JTEMP(M2).E0.JTRCH(M1,1) )  GO=GO+l 

CONTINUE 

IF(JTEMP(M2).EO. JTRCH ( Ml , MLE ) )GO=GO*l 

2  CONTINUE 
IF(GO-l)  3,4,5 

3  MO=MO+l 
MNO(MO)=M] 
GOTO  6 

4  MLF=MLF*1 
ML0RF(MLF)=M1 
GO  TO  6 

5  MB  =  MB-»-l 
MB0TH(MB)=M1 

6  CONTINUE 
1     CONTINUE 

IF  (MB.EO.r)GOTO  7 
DO  10  M5=1,MB 
M6=MBOTH(M5) 
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IF(LENGTH(M6).LT.3)  GOTO  10 

M7=LENGTH(M6) 

IF  (C.NE.O)  GOTO  8 

C=C  +  1 

C0M8(C)=M6 

LENGTH(M6)=-M7 

GOTO  10 

8  CONTINUE 

DO  9  M3=1.C 

M4=C0MB(M3) 

MTEST=JTRCH(M4,1 ) 

DO  13  MB=1,M7 

IF( JTRCH(M6,M8).NE.MTEST>G0  TO  !3 

JD=JD*1 

JDELET( JD)=M6 

LENGTH(M6)=-M7 

GOTO  10 

13  CONTINUE 

9  CONTINUE 
C0MB(C+1)=M6 
LENGTH(M6)=-M7 
C=C  +  1 

ID    CONTINUE 
7     CONTINUE 

IF(MLF.EQ.C  JGOTO  14 

DO  17  Mll=l ,MLF 

M12=ML0PF(M11) 

IF(LENGTH(M12).LT.0)GOTO  17 

M13=LENGTH(M12) 

IF  (C.NE.C)  GOTO  18 

C0MB(C+1)=M12 

LENGTH(M12)=-M13 

C=l 

GOTO  17 
18    CONTINUE 

DO  16  M<?  =  lfC 

M10=C0MB(M9) 

MTEST=JTRCH(M1C,1) 

DO  20  M1A=1,M13 

IF(JTRCH(M12,M14).NE.MTEST)G0T0  2C 

JD=JD+1 

JDELETt JD)=M12 

LENGTH(M12)=-M13 

GO  TO  17 

20  CONTINUE 

16  CONTINUE 
C=C*1 

C0MB<C)=M12 
LENGTH(M12)=-M13 

17  CONTINUE 

14  CONTINUE 
IF(MO.E0.0)GO  TO  21 
CONTINUE 

DO  22  M15=1,M0 
C=C  +  1 
22    C0MB(C)=MN0(M15) 

21  CONTINUE 
RETURN 
END 
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